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SECTION  I 
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INTRODUCTION 

Manned  aircraft  engineering  simulations  are  finding  increased  usage  for  re- 
search and  development  to  reduce  expensive  in-flight  testing.  Industry  and  Government 
organizations  are  developing  multiple  degree-of-freedom  moving  base  simulation 
facilities  which  represent  a significant  capital  investment.  For  efficient  utilization  of 
this  equipment,  methods  must  be  developed  to  plan  and  analyze  simulation  programs 
properly.  Realistic  pilot  models  will  be  one  of  the  foundations  of  these  methods. 

Control  theory  pilot  models  which  accurately  predict  pilot  performance  would 
lead  to  more  efficient  planning,  execution,  and  analysis  of  manned  simulations  of  aero- 
space vehicles.  However,  despite  the  considerable  development  of  control  theoretical 
models  for  the  pilot,  they  have  still  fallen  short  of  adequate  description  of  pilot  behavior 
for  realistic  multi-input,  multi-loop  tasks.  In  particular,  for  these  models  to  be  used  as 
effective  tools  in  simulator  scheduling  and  pre-experimental  analysis,  they  require 
elaboration  to  account  for  the  effects  of  motion  cues,  instrument  and  external  visual 
cues,  and  display  scanning. 

The  Air  Force  Flight  Dynamics  Laboratory  (AFFDL)  has  initiated  a three-phase 
program  to  extend  existing  pilot  modeling  technology  to  more  realistic  pilot  models  for 
application  to  engineering  simulation.  In  support  of  this  program.  Aerospace  Systems, 

Inc.  (ASI)  has  completed  a research  program  on  Phase  1 . Specifically,  this  effort  has 
attempted  to: 

• Assess  past  pilot  models.  I 

• Postulate  an  extended  multivariable  pilot  model  for  a multi-axis, 
multi-input,  multi-output  task  with  motion,  visual  and  scanning 
effects. 


• Demonstrate  the  technical  accuracy  of  the  pilot  model  with 
respect  to  experimental  data  for  a particular  aircraft  and  pilot 
control  task. 

• Document  the  data  used,  the  analysis  methods,  and  the  use  of 
the  pilot  model  postulated. 

Phases  II  and  III,  which  are  not  part  of  the  present  investigation,  will  include 
the  following: 


Phase  II:  Pilot  Model  Verification:  The  Phase  I pilot  model  will 

be  verified  against  a wide  range  of  experimental  data 
for  various  tasks.  This  effort  will  be  developed  through 
a coordinated  AFFDL  and  contractual  program. 

Phase  III:  Simulation  Analysis;  The  verified  pilot  model  will  be 
incorporated  into  a unified  methodology  for  simulation 
planning  and  analysis. 

The  objective  of  the  present  effort  was  to  develop  a multi-variable  pilot 
model  which  can  predict  pilot  behavior  in  realistic  multi-loop,  multi-input,  multi- 
output aircraft  tasks  and  which  also  describes  the  effects  of  motion  cues,  visual  cues, 
and  instrument  scanning.  The  general  pilot  model  developed  during  this  work  should 
improve  simulation  analysis  and  planning,  increase  simulation  efficiency,  and  lower 
simulator  operating  costs.  Furthermore,  aircraft  display  and  control  system  design, 
which  depend  heavily  on  pilot  models,  will  benefit  from  the  results  of  this  effort. 

Our  approach  to  the  development  and  validation  of  the  pilot  model  is  to 
use  data  published  in  the  open  literature.  One  of  the  difficulties  inherent  in  this 
method  is  the  inability  to  obtain  experimental  details  which  might  resolve  any  discrep- 
ancies between  model  predictions  and  experimental  data.  (There  have  been  instances 
where  such  differences  have  lead  to  a discovery  of  anomalies  in  the  experimental  set 
up  rather  than  in  the  model.)  On  the  other  hand,  a model  developed  and  validated 
on  data  from  several  sources  may  prove  to  be  particularly  robust  although  detailed  pre- 
dictions for  any  one  set  of  data  may  not  be  as  precise  as  when  the  data  and  the  model 
come  from  the  same  source. 
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Section  2 of  this  report  presents  j survey  and  assessment  of  past  pilot  models 
for  supporting  manned  simulation  studies.  The  components  of  the  extended  pilot  model 
are  described  in  Section  3,  and  the  development  of  this  model  is  discussed  in  Section  4. 
Section  5 presents  the  results  of  applying  the  model  to  two  previous  experimental  simu- 
lation studies.  Conclusions  and  recommendations  foi  additional  research  are  contained 
in  Section  6.  Details  of  a comprehensive  literature  survey  are  provided  in  Appendix  A. 
The  analytical  formulation  of  the  optimal  control  pilot  model  is  summarized  in  Appendix 
B.  A user's  guide  to  the  computer  program.  Program  PIREP,  which  implements  the 
theoretical  pilot  model,  is  provided  as  a separately-bound  document  for  the  user's  con- 
venience. 


SECTION  II 


SUMMARY  OF  PREVIOUS  RESEARCH 

This  section  presents  an  overview  of  the  previous  literature  concerning  pilot 
models  for  manned  simulation.  We  assume  that  the  reader  is  generally  familiar  with  the 
manual  control  literature  and  with  the  functions  of  the  vestibular  system  for  spatial 
orientation.  (Those  desiring  background  reading  in  manual  control  are  referred  to 
References  1 and  2.  A summary  of  vestibular  functions  is  available  in  Reference  3.) 
Appendix  A contains  a compilation  of  all  pertinent  citations  obtained  during  a compre- 
hensive literature  survey  that  was  conducted  early  in  this  study. 

2.1  GENERAL  CAPABILITIES 

Studies  of  human  operator  behavior  that  have  been  performed  to  date  have 
yielded  a sizable  data  base  from  which  we  may  infer  general,  quantitative  relationships 
between  pilot  model  parameters  and  parameters  of  the  control  environment.  We  can, 
therefore,  use  existing  pilot-vehicle  models  to  predict  system  performance  and  pilot  be- 
havior for  a wide  variety  of  control  systems  and  display  configurations.  Quasi-linear 
models  of  pilot-vehicle  systems  have  been  developed  to  a high  degree.  In  these  models, 
the  pilot  is  represented  by  a linear  response  element  (the  pilot-describing  function)  plus 
a "remnant"  term  (to  account  for  the  portion  of  his  output  that  cannot  be  accounted  for 
by  the  describing  function).  Model  parameters  are  most  readily  interpreied  when  the 
vehicle  dynamics  are  linear  and  when  the  statistics  of  '■he  control  environment  are 
stationary.  The  describing  function  may  then  be  interpreted  as  the  deterministic  linear- 
ized portion  of  the  pilot's  strategy,  with  the  remnant  representing  the  random  non-linear 
and  time-varying  aspects  of  the  pilot's  control.  If  instrument  scanning  is  required,  the 
pilot's  monitoring  strategy  (and  possibly  his  control  strategy)  will  be  time-varying. 
Quasi-linear  pilot-vehicle  models  may  be  usefully  applied  to  situations  of  this  type  as 
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well,  but  measurements  of  pilot  remnant  will  reflect  the  time-varying  nature  of  the 
pilot's  strategy  as  well  as  truly  "random"  behavior. 


Quasi-linear  models  fall  into  two  basic  categories:  the  frequency-domain 
representation  (Reference  1;  and  the  state-variable  (or  "optimal-control”)  model 
(Reference  2).  Both  representations  allow  one  to  predict  pilot  describing  functions  and 
measures  of  overall  system  performance  such  as  mean-squared  system  error.  There  are, 
however,  considerable  differences  in  the  computational  techniques  employed  by  these 
two  types  of  models. 

Both  kinds  of  pilot  models  contain  elements  that  may  be  identified  with 
specific  physiological  or  psychophysical  functions.  The  frequency-domain  models,  for 
example,  often  include  specific  representations  of  the  neuromuscular  system.  On  the 
other  hand,  the  state-variable  model  represents  the  pilot's  estimation  and  control 
strategies  as  distinct  elements.  Frequency-domain  and  state-variable  pilot-vehicle 
models  have  both  been  used  to  explore  the  aspect',  of  pilot  behavior  and  system  per- 
formance described  below. 

• Total  System  Performance.  For  situations  in  which  the  vehicle 
dynamics  can  be  quasi-linearized  and  the  mission  requirements 
represented  in  a suitable  mathematical  form,  analytical  models 
may  be  used  to  predict  overall  system  performance.  Thus,  these 
models  can  be  used  to  determine  the  extent  to  which  mission 
requirements  can  be  met  with  a given  display  configuration.  By 
noting  the  effects  on  system  performance  of  changes  in  one  or 
more  display  parameters,  one  can  (1)  determine  the  information 
requirements  of  the  system,  (2)  determine  improvements  associ- 
ated with  display  quickening  and  predictions,  and  (3)  compare 
the  relative  merits  of  command  information  versus  error  informa- 
tion. Considerable  modeling  success  has  been  obtained  in  this 
regard,  and  the  model  can  now  be  applied  with  reasonable  con- 
fidence to  the  problem  of  predicting  overall  system  performance 
and  how  performance  changes  with  modifications  to  the  display. 

• Pilot  Response  Behavior.  The  input-output  behavior  of  the  pilot 
(i.e.,  the  combined  perceptual,  central-processing,  and  control 
behavior)  is  predictable.  The  response  strategy  is  usually  pre- 
sented in  the  form  of  a describing  function  or  set  of  describing 

1 


-5- 


functions.  Comparison  of  predicted  behavior  with  response 
strategies  measured  in  "conventional"  flight-control  situa- 
tions will  indicate  whether  or  not  an  unusual  response 
strategy  is  demanded  by  the  system  under  investigation. 
Considerable  success  has  been  obtained  in  predicting  pilot- 
describing functions  in  a variety  of  control  situations. 

• The  Pilot's  Estimation  Strategy.  The  pilot's  estimation 
strategy  can  be  predicted  for  situations  in  which  the  pilot 
acts  as  a monitor  only.  A model  of  this  sort  is  quite  valuable 
in  predicting  system  transients  that  occur  when  the  mode  of 
control  switches  from  automatic  to  manual. 

• Pilot  Workload.  Models  have  been  developed  which  appear 
capable  of  predicting  scanning  workload,  i.e.,  the  fraction 
of  attention  that  the  pilot  will  devote  to  the  various  displays 
(References  4,  5).  Models  have  been  applied  to  systems 
containing  physically  separated  displays  as  well  as  those  in 
which  a single  multi-element  display  has  been  used.  In 
addition,  a model  has  been  developed  to  predict  central- 
processing workload,  i.e.,  the  amount  of  attention  that 
must  be  devoted  to  the  task  as  a whole  in  order  for  mission 
requirements  to  be  achieved  (References  5,  6). 

• Instrument-Related  "Noise."  The  state-variable  model  con- 
tains a set  of  parameters  associated  with  perceptual  noise. 

To  some  extent,  these  parameters  may  be  adjusted  to  account 
for  the  effects  of  random  disturbances  inherent  in  the  display 
hardware  (such  as  sensor  and  meter  noise).  They  may  also  be 
used  to  represent  resolution  limitations  of  the  human's  visual 
system  and  random  response  behavior  arising  from  human 
signal/noise  limitations  (References  7 and  8). 

• Handling  Qualities.  Pilot-vehicle  analysis  has  been  applied 
extensively  to  the  study  of  vehicle  handling  characteristics 
(References  9 and  10).  Several  approaches  have  been  ex- 
plored. Although  none  has  achieved  universal  application, 
some  general  observations  can  be  made.  The  pilot's  rating  of 
a vehicle's  handling  qualities  is  influenced  by  the  nature  of 
the  open-loop  response  characteristics  of  the  vehicle,  by  the 
closed-loop  system  performance  that  can  be  achieved  in 
various  task  situations,  and  by  the  nature  of  the  response 
strategy  required  of  the  pilot.  It  has  been  found,  for  example, 
that  pilot  ratings  worsen  as  tracking  error  increases  and  as  the 
lead  that  must  be  generated  by  the  pilot  increases.  Some 
success  has  been  achieved  in  modeling  the  relationship  between 
pilot  rating,  pilot  lead,  and  system  performance.  Since  pilot- 
vehicle  models  are  capable  of  representing  both  the  informational 
characteristics  of  the  display  and  the  perceptual  limitations  of 
the  human,  these  models  should  prove  useful  in  exploring  the 
relationship  between  display  parameters  and  vehicle  handling 
qualities. 
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• Effects  of  System  Failure.  By  investigating  the  effects  of  sudden 
changes  in  the  display  and  system  parameters  on  pilot  behavior 
and  system  performance,  one  can  predict  the  ability  of  the  pilot 
to  recover  from  a failure  condition.  Accordingly,  the  suitability 
of  the  display  configuration  in  a situation  of  this  sort  can  be 
evaluated  (Reference  11). 

• Sensitivity  to  Pilot  Strategy.  The  sensitivity  of  closed-loop 
system  performance  to  changes  in  the  pilot's  response  strategy 
can  be  predicted.  For  example,  one  can  explore  the  effects 
of  a constant  bias  in  one  or  more  of  the  pilot's  response 
parameters,  or  one  can  examine  the  effects  of  an  increase  in 
the  magnitude  of  random  fluctuations  in  these  parameters.  Al- 
though this  aspect  of  the  model's  predictive  capability  cannot 
readily  be  verified  experimentally  (since  the  pilot  will  usually 
attempt  to  select  the  appropriate  response  strategy  with  a mini- 
mum of  random  response  activity),  predictions  of  this  sort  should 
prove  useful  in  the  evaluation  of  system  controllability  and  pilot 
acceptance. 


2.2  MOTION  EFFECTS  ON  THE  PILOT  MODEL 

The  effects  of  vehicle  motion  upon  the  pilot's  control  in  flight  or  in  moving 
base  simulation  has  been  a subject  of  serious  investigation  for  at  least  30  years  (Refer- 
ence 12).  An  accurate  methematical  description  of  the  influence  of  motion  on  pilot 
control  is  a requisite  for  use  of  the  analytical  pilot  model  in  pre-flight  and  pre- 
simulation systems  studies,  and  is  also  necessary  in  the  planning  and  interpretation  of 
aircraft  piloted  simulation  programs,  fixed-base  or  moving-base.  The  increase  in 
recent  years  in  the  capabilities  and  numbers  of  moving  base  simulators  for  research 
purposes  has  intensified  the  effort  to  quantify  the  contribution  of  motion  to  the  pilot 
control  model.  Similarity  and  differences  between  flight  and  fixed-base  simulator 
results  have  been  treated  along  the  following  major  lines: 

• The  additional  lead  generation  produced  by  the  pilot  as  a 
result  of  attitude  motion  sensing  through  his  vestibular 
system,  in  flight  and  moving-base  simulators. 

• The  decrement  in  piloting  performance,  in  flight  or  with 
motion  simulation,  resulting  from  the  direct  effects  of 
sustained  acceleration  and  vibration  upon  the  pilot's 
ability  to  read  instruments  and  perform  fine  motor  control. 
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• The  effects  of  vehicle  moiion  on  a pilot's  perceived  orienta- 
tion in  space,  the  resulting  motion  illusions  associated  with 
this  perception,  and  the  consequent  incidence  of  inappropri- 
ate control  and  aircraft  accidents  associated  with  disorientation. 

In  some  cases  the  thrust  of  the  research  on  motion  cues  has  been  to  use 
physiological  models  of  motion  sensors  to  express  changes  in  the  pilot  control  law  and 
the  resulting  system  performance.  This  has  been  the  approach  taken  in  References  13  - 
15,  all  of  which  use  describing  function  models  for  the  human  operator  and  emphasize 
the  empirical  description  of  motion  effects  on  control. 

In  parallel  with  this  effort,  there  has  been  a major  development  in  mathe- 
matical models  for  the  processing  of  motion  cues  by  the  vestibular  system,  to  produce 
models  for  human  perception  of  attitude  and  position  on  the  basis  of  vestibular  cues 
alone,  and  even  more  recently  on  the  basis  of  combined  visual  and  vestibular  cues. 

The  latest  versions  of  the  vestibular  models  have  been  formulated  in  state  space  nota- 
tion, and  are  based  upon  the  use  of  Kalman  filter  techniques  for  treatmentof  vestibular 
measures  of  pilot  orientation  (Reference  16).  They  are  appropriate  for  incorporation 
in  the  optimum  control  model  of  the  human  operator. 

In  1967,  Young  reviewed  the  results  of  21  previous  investigations  (Reference 
13)  on  comparison  of  pilot  performance  and  pilot  opinion  in  tasks  which  compared 
motion  cues  (simulator  flight)  with  ground-based  fixed  tests.  The  conclusions  were  as 
fol  lows: 

• When  vehicle  dynamics  are  removed  from  the  acceptable  region, 
motion  cues  become  more  helpful  to  the  pilot.  Easily  controlled 
vehicl  es  do  not  require  additional  cues.  Marginally  stable  or 
unstable  vehicles,  however,  may  require  more  lead  compensation 
than  is  comfortably  generated  by  pilots  on  the  basis  of  instruments 
and  consequently  the  rate  sensing  elements  of  the  vestibular 
system  prove  valuable  in  giving  early  indications  of  attitude 
chang  es. 


8 


• Motion  cues  ore  more  helpful  with  unstable  dynamics  than  with 
lightly  damped  oscillatory  vehicles.  As  noted  by  others,  the 
motion  accompanying  high-frequency  vehicle  oscillations  may 
shake  the  pilot  about  in  the  cockpit  making  precise  control  and 
accurate  reading  of  instruments  difficult. 

• Moving-base  simulation  yields  closer  agreement  with  flight 
experience  than  fixed-base  simulation,  especially  in  marginally 
stable  or  difficult  vehicles.  For  most  applications,  simulation 
of  vehicle  angular  motions  is  more  important  than  translation 
simulation.  Stabilization  requires  attitude  control  prior  to 
control  of  position,  thus  emphasizing  rotational  motions.  Normal 
accelerations,  from  centrifuge  simulations,  do  not  appear  as  im- 
portant as  angular  cues  except  under  high  g conditions,  and  the 
spurious  motions  associated  with  closed-loop  centrifuge  simula- 
tion often  cloud  the  pilot  stabilization  picture. 

In  1968  Shirley  (Reference  17)  completed  the  most  comprehensive  investigation 
carried  out  to  that  time  on  the  effects  of  motion  cues  on  the  pilot  model.  Limiting 
himself  to  motion  about  the  roll  axis,  Shirley  compared  the  effects  of  fixed-base  versus 
moving-base  simulation  on  40  different  controlled  elements.  The  range  of  plant  dynamics 
investigated  covers  approximations  to  almost  all  lateral  control  transfer  functions 
encountered  in  flight  vehicles.  Shirley's  results  show  the  differences  in  human 
operator  amplitude  and  phase  angle  over  the  frequency  range  tested  for  roll  motion 
cues  only,  in  combination  with  IFR  visual  cues,  and  finally  for  visual  cues  only. 

In  addition  to  the  gain  and  phase  for  the  pilot  describing  function  as  a function  of 
frequency,  the  remnant  power  spectrum  and  the  measure  of  performance  (mean 
relative  integral  square  error)  is  given.  These  data  form  an  important  source  for 
validation  of  models  which  incorporate  motion  cues.  The  study  shows  how  roll  motion 
cues  permit  the  human  operator  to  increase  his  cross-over  frequency  by  permitting 
higher  gain  tracking  as  a result  of  reduced  phase  lag. 

The  results  of  a number  of  moving-based,  and  fixed-based  simulator 
experiments  in  roll  indicate  that  subjects  can  take  advantage  of  roll  information 
sensed  by  the  semicircular  canals  and  the  otoliths  to  augment  and  substitute  for  visual 


information  in  the  higher  frequency  range  of  vehicle  motion,  but  that  this  information 
cannot  be  substituted  for  visual  information  in  accurate  attitude  or  position  control  in 
a quasi-static  situation  (References  18  - 21). 

In  order  to  distinguish  the  effect  of  linear  acceleration  (sensed  primarily  by 
the  otoliths)  from  the  effects  of  angular  acceleration  (sensed  primarily  by  the  semi- 
circular canals),  Dinsdale  (Reference  22)  investigated  the  relative  effects  of  roll  and 
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yaw  in  manual  control  utilizing  one  of  the  40  cases  of  Shirley's  study  (K/s  ).  By 
examining  the  describing  functions  for  both  otolith  and  canal  stimulation,  he  con- 
cluded that  the  otolith  contributes  an  important  portion  of  the  motion  stabilization 
effects. 

Stapleford,  Peters  and  Alex  (Reference  15)  experimented  with  the  effects  of 
motion  cues  on  the  pilot  describing  function  in  a set  of  roll  control  tasks  with  variable 
lateral  acceleration,  and  were  able  to  identify  separately  the  visual  and  motion  paths 
in  terms  of  effects  on  pilot  describing  function.  Their  conclusions  were  in  general 
agreement  with  those  of  Shirley.  They  found  that  the  visual  path  embodies  higher 
gain  and  less  phase  lead  with  motion  than  without,  presumably  as  a result  of  the  use  of 
vestibular  inputs  to  obtain  roll  velocity  and  acceleration  information.  They  proposed 
models  for  the  utilization  of  vestibular  information  based  on  then  available  mathemati- 
cal descriptions  of  these  sensors:  simplification  of  the  semicircular  canal  path  to  be  a 
first-order  lag  or  even  a pure  time  delay  of  0.1  second,  and  of  the  utricular  path  to  be 
treated  independently  as  a first  order  lag.  Although  these  give  the  correct  direction 
of  effect  on  pilot  describing  function  and  cross-over  frequency,  they  are  inadequate 
for  a pilot  model  intended  to  include  the  effects  of  motion  quantitatively. 

In  a related  study  (Reference  23),  Ringland,  Stapleford  and  Magdaleno 
experimented  with  a multi-axis  VTOL  hovering  task,  with  various  amounts  of  motion 
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involved.  This  study  demonstrated  once  again  the  beneficial  effect  of  roll  motion  and 
the  often  confusing  and  decremental  effect  of  linear  acceleration  in  these  tasks.  The 
results  point  out  the  need  to  consider  the  interaction  between  linear  and  angular 
acceleration  motion  cues  on  the  pilot,  both  in  terms  of  his  perception  of  orientation 
and  his  resulting  control  movements. 

The  performance  effects  of  motion  found  in  the  Stapleford,  et  al.  study 
(Reference  15)  showed  not  only  performance  increase  but  a reduction  in  the  remnant 
component  of  the  mean  square  error  relative  to  fixed  base  values.  Ringland,  et  al. 
(Reference  23)  found,  in  addition  to  the  motion  cue  effects  on  the  describing  function 
discussed  above,  a change  in  the  pilot's  scanning  behavior  relative  to  fixed-base 
conditions.  In  brief,  the  presence  of  motion  cues  permits  the  pilot  to  spend  more  time 
concentrating  on  outer  loop  position  control  (in  terms  of  percentage  of  instrument  scan 
time),  since  higher  frequency  roll  motion  can  be  sensed  by  the  semicircular  canals  and 
otolith s to  help  in  inner  loop  stabilization . These  results  are  consistent  with  those 
found  earlier  by  Benjamin  (Reference  21)  in  fixed-base  and  motion-based  helicopter 
control  studies.  Ringland,  et  al.  (Reference  23)  verified  a large  difference  in  describ- 
ing function  data  between  VMC  (Visual  Meteorological  Conditions)  and  IMC  (Instrument 
Meteorological  Conditions).  This  difference  is  apparent  in  describing  function  mea- 
sures, scan  frequency,  and  pilot  ratings. 

Smaller  differences  in  pilot  function  behavior  between  ground  simulation  and 
motion  (for  this  case,  in-flight  testing)  v/ere  found  by  Newell  (Reference  24)  and 
Newell  and  Smith  (Reference  25).  Newell's  describing  function  data  illustrates  some 
of  the  differences  between  moving-base  simulation  and  flight.  Of  more  practical 
interest  in  Reference  24  is  the  difference  documented  between  IMC  and  VMC  for  both 
flight  and  ground  conditions.  VMC  permit  consistently  larger  amplitude  ratios  for  the 
pilot  and  larger  tracking  errors.  One  must  be  cautious  in  direct  comparison  of  these 
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IMC  and  VMC  in-flight  measurements  of  pilot  describing  functions,  however,  since  they 
were  measured  in  one  case  by  injecting  input  as  a disturbance,  and  in  the  other  by 
injecting  it  as  a command.  In  the  latter  situation,  only  response,  but  not  input  were 
sensed  as  motion  cues.  Qualitatively  similar  results  in  comparison  of  roll  tracking  in 
flight  and  ground  simulation  were  found  by  Smith  (Reference  26).  A summary  of  several 
of  the  earlier  comparisons  of  in-flight  and  simulator  control  for  both  longitudinal  and 
lateral  dynamics  was  given  by  Newell  and  Pietrzak  (Reference  27). 

In  a set  of  experiments  comparing  various  types  of  simulation  and  flight  using 
VTOL  hover  as  the  task,  Sinacori  (Reference  28)  verified  that  the  presence  of  motion 
cues  led  to  the  generation  of  additional  lead  by  the  pilot  and  a marked  decrement  in 
the  tendency  to  overcontrol.  The  smoothness  of  flight  which  is  present  for  cases  involv- 
ing motion,  and  often  absent  in  the  fixed-based  situation,  has  been  noted  by  many 
authors  and  may  be  an  important  resuit  of  the  additional  lead  compensation.  Alterna- 
tively, it  may  be  due  to  a hesitancy  on  the  part  of  the  pilots  to  generate  large 
accelerations  in  a moving-base  simulator  or  in  flight,  as  opposed  to  their  lack  of  such 
inhibitation  in  a fixed-base  simulation  where  the  vehicle  accelerations  are  not  sensed 
directly.  Unfortunately  Sinacori 's  data  includes  information  on  performance  and  stick 
movement  spectral  analysis  but  no  pilot  describing  functions. 

Further  support  for  the  idea  that  accelerations  are  directly  sensed,  and  mini- 
mized, was  provided  by  Rolfe  et  al.  (Reference  29).  By  concentrating  on  the  longi- 
tudinal control  characteristics  they  found  that  when  pitch  motion  was  present  in  a 
moving-base  simulator,  the  control  activity  was  close  to  that  in  flight;  whereas  in  its 
absence,  there  was  a tendency  for  excessive  stick  control  to  be  used,  especially  in 
initiating  and  completing  a turn.  In  addition  to  noting  the  tendency  toward  smoother 
control  in  flight  and  with  motion  cues,  they  note  the  tendency  of  pilots  to  avoid 
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movements  which  would  contribute  to  a negative  g's  situation.  Similar  conclusions  were 
drawn  by  Ringland,  et  al.  (Reference  30).  This  is  furthermore  in  agreement  with  the 
effects  of  motion  on  pilot  control  strategies  found  by  Douvillier,  et  al.,  Fl axman,  and 
Huddleston  (References  31  - 33). 

Bergeron  (Reference  34)  noted,  in  a study  of  pitch  and  roll  axis  control,  that 
in  at  least  some  cases  where  motion  has  relatively  little  effect  on  single  axis  control,  it 
has  a significant  effect  on  both  the  pilot  describing  function  and  the  resulting  system 
error  for  multi-axis  control.  This  presumably  follows  fiom  the  effect  of  motion  in  allow- 
ing the  pilot  to  improve  his  sharing  of  attention  between  two  tasks.  In  an  earlier  study 
of  two-axis  tracking  with  motion,  Bergeron  and  Adams  (Reference  35)  noted  a reduction 
in  the  variability  of  the  measured  human  operator  parameters  when  motion  cues  were 
provided.  This  effect  may  again  be  attributable  to  the  ability  of  the  subject  to  more 
easily  control  tracking  tasks  in  two  separated  axes  when  vestibular  cues  can  be  used  to 
augment  visual  cues  and  alter  scan  patterns. 

Some  substantial  work  of  relevance  to  the  question  of  motion  effects  on 
pilot  models  was  published  during  the  period  of  our  research.  Dillow  et  al.  (Reference 
36)  extended  the  optimal  control  model  for  processing  roll  motion  cues  in  two  ways. 

They  made  roll  acceleration  and  its  derivative  observable,  and  then  placed  weightings 
on  roll  acceleration,  as  well  as  roll  angle  and  rate  in  the  cost  function  to  be  minimized. 
Thresholds  were  assumed  for  roll  and  its  first  three  derivatives.  The  model  extension 
was  tested  against  Shirley's  data  for  a number  of  first,  second  and  third  order  plants, 
under  fixed-base  and  motion-base  conditions.  Using  RMS  roll  error  divided  by  RMS 
input  as  the  sole  performance  measure  (rather  than  pilot  describing  functions).  Reference 
36  demonstrafed  a high  degree  of  agreement  between  analytical  results  and  Shirley's 
experimental  findings. 
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The  most  directly  applicable  new  data  we  uncovered  was  by  Junker  et  al. 
(Reference  37),  and  a companion  paper  by  Repperger  and  Junker  (Reference  38),  These 
papers  treat  tasks  very  similar  to  those  investigated  by  Shirley  and  used  in  our  model 
development,  namely  the  separate  role  played  by  visual  and  motion  cues  in  a roll  track- 
ing task.  The  three  plant  dynamics  used  were  progressively  of  higher  order  and  more 
difficult  to  control,  requiring  more  lead  generation  on  the  part  of  the  pilot.  The  addi- 
tion of  roll  motion  cues  and  the  associated  possibilities  for  lead  generation  resulted  in 
the  largest  amount  of  performance  increment  and  pilot  model  influence  for  the  plants 
in  which  lead  was  most  needed.  The  results  for  several  subjects  were  given  in  terms  of 
RMS  errors,  pilot  transfer  function  frequency  response,  and  values  of  pilot  model  gain, 
and  locations  of  one  zero  and  two  poles  for  each  plant  with  and  without  motion. 


Interestingly,  in  contrast  to  Junker  etal.'s  results,  Shirley  found  that  along 
the  real  axis  of  plant  poles,  for  control  stick  gain  held  constant,  the  motion  cues  were 
of  most  help  for  the  lower  order  plants.  One  must  look  into  the  details  of  the  experi- 
ment to  find  the  basis  for  the  difference.  Shirley's  experiments  were  strictly  compen- 
satory tracking,  in  which  the  command  roll  angle  was  always  zero,  and  the  subject 
rolled  the  simulator  back  to  level  to  overcome  the  disturbances.  He  identified  the 
total  Yp  for  the  case  of  visual  cues  only,  motion  only,  and  vision  plus  motion.  Junker 
et  al.  on  the  other  hand,  used  a roll  tracking  task  in  which  the  pilot  rolled  the  simulator 
to  track  a commanded  roll  angle,  the  error  being  displayed  on  the  instrument.  Thus  for 
the  latter  experiments,  the  visual  tracking  loop  was  strictly  compensatory,  but  the 
motion  feedback  provided  large  angle  actual  response  information.  This  implied 
a pursuit  aspect  to  the  motion  feedback,  and  meant  that  the  actual  angular  rates  and 
positions  sensed  were  much  larger  than  for  the  Shirley  experiments.  The  presence  of 
large  simulator  roll  angles  could  provide  the  additional  stimuli  to  make  the  motion 
sensors  more  effective  in  improving  performance  than  was  the  case  for  the  Shirley 
experiments  with  high  order  systems. 
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Of  particular  value  in  Reference  37  is  the  separate  identification  of  a motion 


cue  transfer  function  acting  on  linear  acceleration  (gravitational  component)  informa- 
tion, and  another  acting  on  angular  acceleration.  The  utricular  information,  based  on 
specific  force,  is  low-pass  filtered  and  utilized  primarily  for  auasi-static  position 
sensing.  (Junker  et  al.  found  no  lead  generation  in  the  linear  acceleration  pathway.) 
On  the  other  hand,  the  semicircular  canals,  detecting  angular  acceleration,  are  used 
extensively  for  roll  angle  lead  generation  when  they  are  stimulated.  (Junker  et  al. 
find  that  the  angular  acceleration  path  transfer  function  changes  from  low  pass  to  high 
pass  with  motion  cues  in  the  higher  order  plants.) 

In  addition  to  the  pilot  model  results.  Junker  and  Replogle  (Reference  39) 
also  demonstrate  the  increased  effectiveness  of  roll  motion  cues  in  reducing  task  learn- 
ing time  and  improving  tracking  performance  for  tasks  having  large  roll  motion  and 

2 

requiring  lead,  using  a plant  of  K/s  . 

2.3  SIMULATOR  REQUIREMENTS  FOR  MOTION  CUES 

In  addition  to  the  research  directed  precisely  at  the  question  of  motion  cue 
effects  on  pilot  performance,  there  has  been  considerable  interest  in  determining  the 
adequacy  and  importance  of  simulator  motion  and  the  development  of  optimum  wash- 
out circuits  for  a given  simulator  configuration.  A review  of  the  factors  in  simulator 
design  which  relate  to  motion  and  visual  cues  is  given  by  Matheny  et  al.  (Reference 
40).  They  generaUy  support  the  view  of  the  pilot  motion  cue  utilization  discussed 
above.  For  attitudinal  changes,  motions  will  be  sensed  first  through  the  vestibular 
senses.  The  visual  senses  are  generally  relied  upon  to  corroborate  or  substantiate  the 
information  received  by  the  motion  senses,  although  in  some  systems  the  rote  of  move- 
ment is  such  that  the  visual  senses  are  able  to  detect  the  fact  of  movement  before  the 
motion  senses. 
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In  an  experimental  study  of  factors  in  flight  simulation  aimed  at  the  required 
fidelity  of  motion  cues,  Klier  and  Gage  (Reference  41)  studied  a simplified  air-to-air 
combat  tracking  task  in  a moving  base  simulator  with  pitch,  roll  and  heave  degrees  of 
freedom.  Using  various  cut-off  frequencies  for  simulator  motion,  they  found  the  best 
performance  in  terms  of  time-on-target  and  RMS  error,  with  a 0 to  2.5  Hz  bandwidth. 
Performance  dropped  as  the  bandwidth  was  increased  or  decreased  from  this  level.  The 
results  obtained  by  Klier  and  Gage  support  the  idea  that  motion  cues  lead  to  vestibular 
generation  of  lead  with  good  tracking  at  frequencies  below  2 Hz,  but  that  at  higher 
frequencies  the  presence  of  vibration,  which  cannot  be  damped  out  by  the  pilot, 
serves  only  to  degrade  vision  and  interfere  with  his  ability  to  track.  In  another  study 
using  two-axis  movement  control  in  a moving-base  and  a similar  flight  base  simulation 
for  pilot  training,  Jacobs  et  al.  (Reference  42)  verified  the  importance  of  simulator 
motion  in  pilot  steering  performance,  and  found  an  important  interaction  between 
simulator  motion  and  preference  for  specific  display  formats. 


In  summarizing  the  state  of  art  in  this  area,  Staples  (Reference  43)  pointed 
out  the  importance  for  proper  vestibular  modeling  and  its  use  in  simulator  wash-out 
design.  In  his  discussion  of  Staples'  paper,  Bray  (Reference  44)  noted  that  large  dif- 
ferences in  the  pilot's  estimation  of  spatial  orientation  cun  result  from  their  motion- 


visual  interaction,  and  that  lack  of  motion  cues  causes  a much  large:  decrement  in 
performance  for  those  subjects  classified  (on  a rod  and  frame  tert)  as  •'frame  oriented" 
than  for  those  who  were  not  frame-oriented.  This  discussion  further  emphcnizes  the 
importance  of  systematic  investigation  of  visual-vestibular  interaction  and  spatial 
orientation,  both  in  terms  of  how  the  pilot  feels,  and  the  way  he  will  react  in  a situa- 
tion tending  to  produce  false  motion  illusions.  Some  of  the  same  thoughts  are  expressed 
by  Lew  and  Dyda  (Reference  45)  in  a discussion  of  the  pros  and  cons  of  ground-teased 
simulation  and  comparison  with  flight  for  landing  training. 
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Sinacori  (Reference  46)  discusses  the  need  for  realistic,  distortion-free  visual 
displays  in  a VFR  moving-base  simulation.  He  also  lays  out  the  goals  of  their  motion 
drive  system  in  pitch,  roll  and  yaw  to  provide  the  following: 

1 . Provide  semicircular  canal  stimulation  at  frequencies  where 

the  stimulus  corresponds  nearly  to  perceived  angular  velocity, 

i.e.,  0.2  to  10  rad/ sec . 

2.  Provide  utricular  stimulation  at  low  frequencies  only,  i.e., 
up  to  0.25  rad/sec. 

3.  Minimize  false  stimulation  of  the  vestibular  organs  and  main- 
tain low  frequency  postural  reflexes. 

2.4  DIRECT  EFFECTS  OF  SIMULATOR  MOTION  ON  PILOT  PERFORMANCE 
INDEPENDENT  OF  MOTION  CUES  UTILIZED  IN  FEEDBACK 

It  is  well  known  that,  in  addition  to  the  beneficial  effect  of  simulator  or 
flight  motion  on  improving  performance,  aiding  pilot  rating,  and  leading  to  smoother 
controls,  the  presence  of  vehicle  motion  in  itself  can  lead  to  decrements  in  performance, 
especially  when  the  motion  is  sustained  high  acceleration  or  of  high  frequency  vibration. 
Some  of  the  high  acceleration  effects  on  the  pilot  describing  function  were  reviewed  by 
Sadoff  and  Dolkas  (Reference  47).  Recent  investigations  of  the  effect  of  vehicle  vibra- 
tion on  pilot  describing  functions  were  reviewed  at  the  Tenth  Annual  Manual  Control 
Conference  ^References  48  and  49), 

Possibly  less  well  known  as  a direct  effect  of  motion  on  pilot  performance 
are  experiments  performed  by  Guedry  and  colleagues  at  Pensacola.  Benson  and  Guedry 
(Reference  50)  demonstrate  conclusively  the  effect  of  rotary  vehicle  motion  in  decreas- 
ing the  performance  on  an  unrelated  tracking  task  as  mediated  through  the  oscillatory 
movement  of  the  eyes  called  nystagmus.  Larger  performance  decrements  are  associated 
with  larger  values  of  nystagmus  in  both  the  pitch  and  yaw  axes.  Gilson,  Benson,  and 
Guedry  (Reference  51)  in  the  same  year  show  that  as  instrument  illumination  is  increased 


him.,..,  m 

to  reduce  the  level  of  nystagmus  during  vehicle  rotation,  the  interference  with  tracking 
is  diminished.  Both  of  these  contributions  lead  to  the  conclusion  that  an  important 
direct  effect  of  motion  on  pilot  control  is  mediated  through  the  development  of  nystagmus 
and  its  resulting  impairment  of  visual  acuity.  This  can  be  represented  in  the  human 
operator  model  as  an  increase  in  observation  noise. 


- 18  - 


SECTION  III 


COMPONENTS  OF  THE  PILOT  MODEL 

This  section  discusses  the  principal  components  of  the  pilot  model  formulated 
for  manned  simulation  applications.  The  optimal  control  model  of  the  human  operator 
is  selected  as  the  foundation  for  this  extended  model . Workload  and  scanning  effects  are 
described,  followed  by  a discussion  of  motion  cues  and  the  vestibular  sensors.  The  dif- 
ferences between  visual  and  instrument  flight  conditions  are  then  explained.  Finally,  a 
summary  of  the  extended  pilot  model  is  presented  for  convenient  reference. 

3.1  INTRODUCTION 

As  discussed  in  Section  2,  there  are  two  basic  types  of  pilot  models:  describ- 
ing function  models  and  state  space  models.  The  first  type,  which  is  formulated  in  the 
frequency  domain,  originated  with  classical  control  theory.  It  has  been  developed  to 
control  various  single-  and  multi-loop  system  configurations.  The  second  type  of  pilot 
model,  which  is  formulated  in  the  time  domain,  was  developed  from  modern  control  and 
estimation  theory.  It,  too,  is  applicable  to  a variety  of  system  configurations.  These 
approaches  to  pilot  modeling  have  been  effectively  applied  in  the  fields  of  aircraft  fly- 
ing qualities,  pilot-vehicle  control  system  integration,  and  display  synthesis. 

The  optimal  control  model  of  the  human  operator  has  been  selected  in  this 
research  effort  for  the  following  reasons: 

• The  model  more  easily  handles  multi-input  and  multi -output 
control  tasks. 

• The  model  provides  an  empirically  verified  measure  of  workload 
which  is  related  to  attention  rather  than  eye  movement  behavior. 

As  such,  it  will  also  work  at  the  level  of  displayed  information  or 
the  lavel  of  detailed  display  format  for  both  integrated  and 
separate  instrument  displays. 
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• The  optimal  control  model  of  the  human  operator  is  more  adaptable 
to  calculating  the  time  varying  statistical  behavior  (over  the  en- 
semble of  trajectories)  due  to  nonstationary  statistics  arising  from 
the  geometric  and  kinematic  relationships  of  the  guidance  and 
navigation  systems. 

In  the  remainder  of  this  section  we  briefly  outline  the  baseline  optimal 
control  model  of  the  pilot,  discuss  the  extensions  which  incorporate  scanning  and  atten- 
tion allocation,  motion  cues,  and  VMC/IMC  (Visual  Meteorological  Conditions/ 
Instrument  Meteorological  Conditions)  cues.  Lastly,  we  summarize  the  parameters  of 
the  models,  guidelines  on  how  to  choose  the  parameters,  and  what  effects  are  repre- 
sented by  the  parameters. 

3.2  THE  OPTIMAL  CONTROL  MODEL  OF  THE  PILOT 

This  section  briefly  reviews  the  important  elements  of  the  optimal  control 
model  of  the  pilot  in  describing  human  performance  for  fixed-base  IMC  (instrument 
meteorological  conditions).  These  conditions  are  the  ones  for  which  the  model  has  been 
well  validated  in  the  past.  A more  complete  analytical  description,  including  the 
relevant  equations  of  the  model,  is  contained  in  Appendix  B. 


A block  diagram  of  the  model  is  presented  in  Figure  1,  which  shows  the 
information  flow  and  the  interaction  of  the  pilot  with  the  controlled  element.  The  dis- 
played information  y(t)  passes  through  the  threshold  element  (either  a psychophysical 
threshold  or  an  indifference  threshold).  Perceptual  observation  noise  v (t)  is  added  to 
this  quantity,  which  then  passes  through  the  perceptual  time  delay  t.  The  perceived 
displayed  quantities  y^(t)  are  the  inputs  to  the  core  of  the  model,  which  consists  of  a 
Kalman  filter  to  account  for  the  controlled  element  dynamics  and  the  observation  noise, 
and  a predictor  to  compensate  for  time  delay.  The  current  state  estimate  x(t)  is  com- 
bined linearly  to  provide  a commanded  control  signal.  To  this  control  signal  is  added 
a white  motor  noise  component  v^t)  with  the  sum  passing  through  a neuromuscular/ 


20 


Figure  1 . Optimal  Control  Model  of  Human  Response. 


manipulator  first-order  lag.  For  a single  displayed  quantity,  the  covariance  of  the 
observation  noise  depends  on  several  important  parameters  and  has  the  form 


V 

*i 


(1) 


where  in  the  above  expression  V is  the  covariance  of  the  white  observation  noise; 

• th 

a.  is  the  standard  deviation  of  the  i displayed  quantify;  N is  the  describing  function 
for  the  dead  zone  element  of  half  width  o.;  and  p.  is  the  multiplicative  noise/signal 
ratio. 


For  fixed-base  IMC  tracking,  it  has  been  found  that  a noise/signal  ratio  of 
-20  dB  (or  0.01  n)  provides  a good  match  of  the  model  with  the  experimental  data.  It 
has  also  beer  ound  that  for  each  displayed  element,  there  should  be  two  inputs  to  the 
model  corresponding  to  the  position  of  the  display  element  and  its  rate.  The  above 
relationship  on  noise  covariance  should  be  used  for  each  of  these  Kalman  filter  inputs. 

One  element  of  the  model  which  is  not  shown  in  the  block  diagram  is  the 
quadratic  functional  which  is  minimized  by  the  choice  of  the  control  gains,  L,  shown 
in  Figure  1.  In  its  most  general  form  the  functional  is  given  by 

T 

J = E | lim  f [y'Q  y + x‘Q  x + u'Q  u + u'Q  u]  dt 
IT-*  oo  J Y * r u 
0 

where  Q , Q^,  Q , and  are  symmetric  weighting  matrices  for  the  display  variables, 
state  variables,  control  rate,  and  control  displacement,  respectively. 

The  specific  choice  of  these  weightings  will  be  discussed  later,  but  in  prac- 
tice it  has  usually  been  found  sufficient  to  weight  the  display  deviations  and  control 
rate  only,  with  the  weights  on  control  rate  chosen  to  provide  a neuromuscular  time 
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constant  (T^  in  Figure  1)  consistent  with  the  manipulator  (with  a value  of  - 
0.1  - 0.2  second  being  typical). 

3.3  WORKLOAD  AND  SCANNING  MODELS 

Estimates  of  pilot  workload  developed  to  date  have  been  of  two  major  types: 
visual  scanning  workload  (Reference  4)  or  task  interference  workload  (Reference  6). 

Both  have  been  well  documented  and  verified  in  experimental  situations.  There  are 
several  reasons  why  we  elect  the  task  interference  workload  model  in  this  study. 

• Its  formulation  is  consistent  with  the  optimal  control  model' 
of  the  pilot 

• It  will  predict  time  variations  and  workload  level  in  non- 
stationary conditions 

• It  may  be  used  either  at  the  information  requirements  level 
or  the  display  format  level  of  detail 

• It  contains  flexibility  to  accommodate  side  task  workloads. 

We  feel  this  last  item  is  of  particular  importance  in  predicting  pilot  workload  models  for 
the  purpose  of  simulator  resource  allocation.  Many  of  the  pilot  model  applications  to  date 
have  concentrated  solely  on  the  continuous  control  tasks.  However,  it  is  known  that  there 
are  many  side  tasks  of  a discrete  nature  (e.g.,  communications,  frequency  changes,  equip- 
ment adjustments,  etc.)  which  aie  not  usually  accounted  for  in  most  modeling  efforts, 
but  which  can  be  of  great  importance.  The  visual  scanning  method  of  measuring 
pilot  workload  appears  not  to  have  the  advantages  listed  above  because  it  is  based  on 
time  averaged  (stationary)  statistics, and  the  scanning  behavior  is  too  detailed  a metric 
when  working  at  the  information  requirements  level  of  detail  (References  51-54). 

We  let  the  total  fraction  of  attention  (nominally  1.0)  be  composed  of  the 
additive  terms 
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where 


f,  = f + f + f „ 

tot  c o other 


f = r f 

c . c. 


In  the  above  equation,  is  that  portion  of  the  total  attention  devoted  to  the  control 

th 

task  (control  task  workload);  f is  the  fraction  of  attention  allocated  to  the  i display 

ci 

variable  for  control;  fQ  is  the  fraction  of  attention  lost  in  switching  from  one  display  to 
another,  or  from  control  task  to  side  task;  and  f ^er  's  the  fraction  of  attention  allo- 
cated to  the  other  side  tasks  demanding  attention  from  the  pilot. 

The  effect  of  any  given  allocation  of  control  task  workload  among  the  dis- 
play variables  (i.e.,  specifying  the  individual  f ) has  an  effect  on  performance  as 

i 

follows.  The  task  interference  measure  of  workload  indicates  that  for  a modification 
of  the  covariance  of  observation  noise  of  each  input,  a Kalman  filter  in  the  pilot  is 
necessary  (References  6,  55-57): 


V = - 
y;  f 


y;  _ 


fc.  LN-ai'  a\) 

i 


0 fh 

where  VU  is, the  power  spectral  density  of  the  observation  noise  of  the  i display 

yi 

variable  when  attention  is  limited  to  that  display  variable  alone  (Equation  (1)),  and 

t fl 

f is  the  fraction  of  attention  allocated  to  the  i display  variable.  This  workload 
c • 

I 

measure  has  the  correct  interpretation  at  the  limiting  values  of  full  attention  (fc  = 1) 

and  no  attention  (f  = 0). 

c. 

Using  this  model,  we  may  determine  the  attentional  workload  requirements  as 

a function  of  system  performance  as  follows.  Assume  a given  value  of  fraction  of  attention 

on  the  control  task  workload  t . Then  determine  the  f within  the  constraint  of  Equation  (4) 

c c. 
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to  minimize  the  quadratic  performance  index  (see  Appendix  8).  There  arc  the 

optimal  allocations  of  attention  between  display  variables  for  a specified  fraction  of 

attention  allocated  to  the  control  task.  From  the  RMS  performance  calculations,  we 

also  know  the  system  performance  J.  If  this  process  is  performed  for  several  different 

values  of  f , then  we  may  draw  the  curve  o'  control -task  workload  required  to  accom- 
c 

plish  the  given  level  of  system  performance  (J)  as  shown  in  Figure  2.  The  solid  curve 
represents  the  control  task  workload,  and  the  dashed  curve  represents  the  effect  of  add- 
ing the  fraction  of  attention  for  transitions  and  the  fraction  of  attention  required  for 
other  tasks.  This  dashed  line  represents  the  total  required  fraction  of  attention  for  the 
pilot  and  represents  the  performance  workload  tradeoffs.  Note  that  the  pilot  is  work- 
ing at  full  capacity  to  maintain  the  cost  functior.  at  or  below  Jj . 


Figure  2.  Required  Attentional  Workload  vs  System  Performance. 
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MOTION  CUES:  STATE  SPACE  FORMULATION  OF  THE  VESTIBULAR 
SENSORS 


In  this  section  we  formulate  the  equations  to  be  used  in  the  optimal  control 
model  of  the  pilot  to  represent  the  vestibular  sensors.  It  consists  of  two  sections:  the 
development  of  the  semicircular  canal  equations,  and  the  development  of  the  otolith 
equations. 

3.4.1  SEMICIRCULAR  CANALS 

In  Equation  (3.9)  of  Ormsby's  thesis  (Reference  16),  the  Laplace  transform 
equation  for  the  afferent  firing  rate  of  a semicircular  canal  is  given  by 


AFR(s) 


300s  (.02s  + 1 ) 


(impulses/second)  ( 1 8s  + 1 ) (30s  + 1 ) 


u>(s ) + v (s  ) 


where  the  input  angular  acceleration  and  additive  noise  satisfying 

2 

u>(s)  = angular  acceleration  (deg/sec  ) 


vscc(0  Vscc(t)  = (.223)z  6(t - r) 


In  the  above  equation  we  have  used  a rate  sensitivity  parameter  of  .02  seconds 
(rather  than  the  .01  seconds  of  Ormsby)  because  this  seems  to  be  a better  fit  to  the 
data;  we  have  also  neglected  the  very  short  time  constant  (.005  seconds).  This 


equation  may  be  rewritten 


AFR(s)  = .01111  1 + 


49.911s  - .001852 
s2  + .08888  +.001852 


w(s ) + v. 


Ormsby's  equations  are  calibrated  in  terms  of  impulses  per  second,  and  we 
prefer  to  scale  the  output  in  units  of  angular  acceleration.  To  do  this,  we  use  the 
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fact  that  near  the  psychophysical  threshold  (taken  to  be  0.19°/sec  ),  the  ratio  of 

2 

signal  to  noise  variance  is  approximately  (1 .36)  to  provide  the  probability  of 

2 

correct  detection  of  .75.  Thus,  in  the  scaled  equation  the  noise  variance,  , 

9 

in  units  of  (degrees/second  *)  is  given  by 


(.  19°/sec2)2  l1’36 


(8) 


Using  the  fact  that  the  output  will  be  scaled  by  the  factor  K,  then  we  have  the 
relation 


where 


2 - *2  2 

a . — K a 

v v 


(9) 


see 


K = 


0. !9°/sec2 
(1.36)  (.223)  ips 


Thus  we  are  finally  led  to  the  following  state  equations  and  measurement  equations 
for  the  semicircular  canals 


Xsccl  Xscc2 


x « = -.001852  x , - .08888  x , 
scc2  sccl  see? 


♦ w 


(10) 


yscc  ^e9/ sec3)  " -1.289  X 10  ® . 

* 0.00696  w ♦ v 


*»cc2 


(H) 
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3.4.2  OTOLITHS 


Ijie  state  equations  and  measurement  equations  for  the  otolith  are  taken  from 
Ormsby's  thesis  (Reference  16)  and  are  aiven  by 


x -.2  x + .2  f 

oto  oto 


y - -46 .lx  + 9 1 . 1 f + v , 
'oto  oto  oto 


where  f is  the  specific  force.  We  again  are  faced  with  the  output  equation  in  units 
of  impulses  per  second.  Ormsby  uses  the  scale  factor  such  that  45  impulses  per  second 
corresponds  to  a steady  state  stimulation  of  lg  . Thus  we  may  rescale  the  output 
equation  to  be 


y (in  g's)  = -1 .024  x + 2.024  f + v 
'oto  a oto  oto 


voto(0  vQto(r)  - (.003267  «,)*  sec  6(t  -r) 


3.5  VMC/IMC  CUES 

One  of  the  goals  (of  this  research  effort  is  to  account  for  outside-the- 
window  visual  cues  and  their  effects  on  pilot  behavior  in  comparison  to  the  normal 
instrument  (IMC)  cues.  This  study  was  limited  to  the  effect  on  attitude  control  per- 
formance, and  did  not  consider  the  more  complicated  position  control  task  such  as 
occurs  during  the  final  phase  of  landing  (i.e.,  near  flare  height  and  below).  The 
effect  of  full  field  visual  cues  has  been  measured  in  flight  and  in  ground-based  simula- 
tion using  a wide  angle  contact  analog  display  (Reference  25).  The  effects  of  the  wide 
field  display  manifested  themselves  in  higher  pilot  describing  function  gains  when 
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compared  to  the  IMC  tracking  tasks,  both  with  and  without  motion.  These  effects  might 
be  accounted  for  in  the  model  by  decreasing  the  noise/signal  ratio  on  attitude  measurement, 
especially  attitude  rate;  by  decreasing  the  threshold  upon  the  corresponding  measurements; 
and  by  decreasing  the  weight  on  the  attitude  rate  in  the  quadratic  cost  functional.  All 
of  the  e would  have  the  effect  of  increasing  pilot  gain  (although  to  varying  degrees)  and 
it  is  one  of  the  purposes  of  our  investigation  to  determine  the  most  appropriate  way  to 
model  the  VMC  cues. 

3.6  SUMMARY  OF  THE  PILOT  MODEL 

In  this  section  we  review  the  parameters  of  the  optimum  control  model  and 
briefly  describe  how  one  chooses  the  parameters  to  represent  the  situation  under  study. 

The  outline  of  this  discussion  will  closely  follow  the  description  of  the  program  input  con- 
tained in  Volume  2 of  the  report. 

3.6.1  CONTROLLED  ELEMENT  DYNAMICS:  MATRICES  A,  B,  E 

These  matrices  for  the  system  response  (See  Appendix  B)  describe  the  dynami- 
cal characteristics  of  the  element  being  controlled  and  any  shaping  filters  required  to 
describe  the  spectral  characteristics  of  the  disturbance  input.  In  the  use  of  the  model 
to  describe  motion  cues,  these  matrices  should  also  reflect  the  vestibular  dynamics  dis- 
cussed in  Section  3.4. 

3.6.2  COST  FUNCTIONAL  WEIGHTINGS 

The  cost  functional  weighting  matrices  contain  only  diagonal  elements  in  the 
usual  application  of  the  optimal  control  model  of  the  human  operator.  In  addition,  it 
has  generally  been  found  sufficient  to  ignore  weightings  on  the  state  and  on  the 
control  displacement  Q^.  The  weightings  on  the  displayed  quantities  and  a control 
rate  weight  Qr  are  discussed  next 
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We  have  generally  obtained  excellent  results  from  the  model  by  weighting 
the  position  of  the  displayed  variables  inversely  proportional  to  the  square  of  the  maxi- 
mum allowable  displacement.  When  dealing  with  a specific  display  format  layout,  this 
generally  corresponds  to  full  scale  on  the  display.  When  dealing  with  the  model  at  a 
more  general  level,  this  would  correspond  to  a maximum  display  value  appropriate  to 
specified  mission  requirements.  Recall  that  the  optimum  control  model  of  the  pilot  re- 
quires two  inputs  for  each  displayed  quantity:  its  position  and  its  rate.  While  the  dis- 
cussion above  describes  how  to  choose  the  weight  on  the  position  of  the  display  element, 
there  has  previously  been  little  guidance  available  for  selecting  the  weight  for  the  rate 
of  a displayed  quantity. 

In  many  applications  the  maximum  value  of  a position  or  attitude  variable  is 

easy  to  determine  from  the  mission  requirements,  hence  its  weighting  coefficient  in  the 

-2 

quadratic  cost  functional  is  relatively  straightforward,  i.e.,  the  coefficient  is  (ymax) 
There  is  no  corresponding  interpretation  on  the  rate  of  change  of  that  variable  in  many 
instances,  so  it  is  subject  to  some  uncertainty.  Figure  3 shows  the  change  in  RMS  error 

O 

predicted  by  the  model  as  a function  of  (y  /y  ) for  several  of  the  plants  treated 

inox  max 

2 

by  Shirley  (Reference  17).  The  abscissa  has  units  of  seconds  , and  its  interpretation 
led  to  our  final  choice  of  the  weighting  coefficients.  As  shown  in  Figure  3,  this  can 
lead  to  significant  changes  in  the  performance  predictions  of  the  model.  To  place  a 
realistic  value  on  the  weighting  of  the  rate  term,  consider  that  if  the  minimization 
process  is  doing  a reasonable  job,  we  would  expect 
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Figure  3.  Comparison  of  Model  Predictions  with  Change  in 
Weighting  Coefficients  - Visual  Only. 


If  this  ratio  remains  true  for  all  levels  of  y and  y,  and  if  the  system  exhibits  stable 
behavior,  then  Equation  (16)  can  be  approximately  written 

TRy  + y»l 


where  the  equivalent  time  constant  of  the  response  is  given  by 


Thus,  given  a value  of  ymax/  and  the  system  requirements  for  the  response  time,  one 

can  choose  a reasonable  value  of  y for  the  quadratic  weightings.  For  example,  in 

max 

attitude  control  tasks  where  a typical  response  time  would  be  on  the  order  of  2-3  seconds, 

the  t jivalent  time  constant  (td)  should  be  on  the  order  of  1 so  that  y is  approximate- 
1 K max 

ly  equal  to  y .In  trajectory  control  problems,  the  response  time  will  be  considerably 
max 

longer,  hence  y is  much  less  than  y . 
a ' 'max  'max 

The  other  weighting  coefficients  in  the  cost  functional,  those  for  control  rate, 
are  generally  chosen  (for  a well-designed,  lightweight,  force-type  manipulator)  to  pro- 
vide a neuromuscular  time  constant  of  0.1  second.  There  are  many  situations  where 
this  may  not  be  the  correct  value  to  use,  however.  For  example,  if  the  manipulator  is 
a large  displacement,  high  force  gradient  device  such  as  elevator  control  on  a large 
transport  aircraft,  then  neuromuscular  lags  considerably  larger  than  0.1  second  should 
be  expected  to  be  more  representative,  e.g.,  0.2  to  0.5  seconds.  In  addition,  there 
is  an  absolute  limit  on  the  velocity  of  hand  movement  that  can  be  expected,  with 
values  of  2 - 5 inches  per  second  being  representative  for  a well-designed  manipulator. 
Thus,  the  control  rate  weighting  factor  is  chosen  to  provide  the  specified  neuromuscular 
time  constant;  if  this  value  corresponds  to  control  rates  larger  than  2-5  inches  per  second, 
then  this  latter  value  should  be  considered. 

3.6.3  CONTROLLER  TIME  DELAY 

This  parameter  represents  the  pure  transport  time  delay  of  the  human  operator, 
but  also  can  be  used  to  model  the  transport  delays  in  the  controlled  element.  The 
typical  value  for  human  operator  reaction  times  is  0.2  second. 


J 
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3.6.4  PROCESSOR  DRIVING  NOISE 


This  is  the  covariance  of  the  white  process  noise  driving  the  state  equations. 
It  is  chosen  to  provide  the  appropriate  RMS  disturbance  value. 

3.6.5  ADDITIVE  MOTOR  NOISE 

This  is  the  covariance  of  the  white  noise  component  added  to  the  human 
operator  output  (See  Figure  1).  It  is  typically  chosen  to  provide  a noise/signal  ratio 
referenced  to  the  control  output  of  -25  dB  (0.003^),  but  there  are  situations  where 
this  may  lead  to  unrealistically  small  noise  levels.  In  addition,  poorly  designed 
manipulators  with  much  stiction  and/or  the  effects  of  motion  may  increase  this  above 
the  -25  dB  level . 

3.6.6  OBSERVATION  THRESHOLDS 

These  parameters  of  the  model  generally  are  chosen  to  represent  indifference 
thresholds  and  are  strongly  dependent  on  the  format  of  the  display.  Typically,  a 
reasonable  value  for  these  thresholds  is  the  minimum  resolution  marking  on  the  display, 
but  in  other  instances  the  more  appropriate  value  may  be  (say)  10  percent  of  the  full- 
scale  value.  A pilot  opinion  is  useful  in  this  context  also,  since  one  can  assess  the 
indifference  thresholds  by  asking  experienced  pilots  what  type  of  a static  deviation 
they  are  willing  to  tolerate. 

The  threshold  for  the  display  element  rate  input  can  be  chosen  in  a manner 
analogous  to  the  weighting  of  the  same  quantity,  i.e.,  for  an  equivalent  time  constant 
of  response  we  would  expect  that 


a 
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where  a and  ».  are  the  threshold  halfwiu  ''$  on  the  displayed  quantity  y and  its  rate  y. 

y y 

That  is,  the  ratio  of  the  threshold  of  display  position  to  display  rate  is  the  same  as  the 
ratio  of  maximum  displayed  quantity  to  its  maximum  rate. 

These  thresholds  enter  the  calculation  of  the  observation  noise  in  the  pilot 
model  through  the  equation  (See  Equation  (1)  ) 

P. 

V = -L_ 

y,  f 

Cj 

where  in  the  above  expression,  N is  the  describing  function  of  the  dead  zone  element 
shown  in  Figure  1,  and  ar.  is  the  threshold  value  described  above. 

3.6.7  OBSERVATION  NOISE  RATIOS 

The  observation  noise  ratios  are  the  parameters  P . in  Equation  (20)  above, 
and  for  visual  displays  it  has  been  determined  that  a value  of  0.01  n (-20  dB)  generally 
provides  a good  fit  with  experimental  data.  However,  for  motion  cue  inputs  (and  per- 
haps VMC  cue  inputs)  other  values  of  noise/signal  ratios  may  be  appropriate.  This  ratio 
will  be  one  of  the  parameters  which  will  be  varied  in  the  examination  of  motion  cue 
information. 

3.6.8  ATTENTION  ALLOCATION 

The  attention  allocations  are  the  fractions  of  attention  devoted  to  each  dis- 
play element  subject  to  the  constraint  of  Equation  (4).  The  fraction  of  control  atten- 
tion devoted  to  the  display  enters  into  the  computation  of  the  noise  covariances  as 
indicated  in  Equation  (20).  It  is  through  this  effect  that  we  expect  to  see  differences 
between  scanning  behavior  in  motion  and  no-motion  conditions.  The  vestibular  inputs 
will  be  assigned  f =1.0,  and  the  fraction  of  attention  for  the  visual  displays  will  be 

Ci 

chosen  to  further  minimize  the  performance  index  as  indicated  in  Section  3.3. 
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N(\,  / a: ) 


(20) 


In  its  most  general  form  the  mod  I is  truly  a multi-input,  multi-axis  representa- 
tion. However,  when  dynamical  decoupling  exists  between  axes,  the  attention  alloca- 
tion can  be  optimized  on  each  axis  separately  for  a few  values  of  f . Then  a final 
optimization  can  be  easily  performed  to  determine  the  overall  allocation  of  attention 
among  axes.  This  procedure  provides  considerable  saving  in  computational  effort  com- 
pared to  an  attention  allocation  among  all  displays  simultaneously. 

3.6.9  SUMMARY 

Table  1 summarizes  the  correspondence  of  the  real  world  components  and 
the  model  components  as  they  exist  at  this  point  of  the  discussion.  Note  that  the 
VMC/IMC  cue  may  be  represented  by  functional  weightings,  observation  noise  ratios, 
or  thresholds.  This  ambiguity  will  be  clarified  in  later  sections  of  the  report  when  we 
determine  the  relative  sensitivity  of  model  predictions  to  these  parameters  for  VMC/ 

IMC  cues. 
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Table  1.  Correspondence  of  Real  World  and  Model  Components. 


Real  World  Component 

Model  Component 

Airframe  Dynamics 

Augmentation  System 

Mission  Phase 

Vestibular  Afferent  Dynamics 
Simulator  Dynamics 

Control led-element  Dynamics 

x = Ax  + Bu 

Displays 

(Visual  and  Vestibular) 

Measurement  Vector 

y = Cx  + Du 

Mission  Requirements 

VMC/IMC  Cues 

Neuromuscular/Manipulator 

Dynamics 

Cost  Functional  Weightings 

J = Ej  f[yTQyy+uTQrUj  dtj 

Vestibular  Cues 

VMC/IMC  Cues 

Observation  Noise/Signal  Ratio 
Display  Thresholds 

Disturbance  Environment 

(Turbulence,  Shear, 

Guidance  Errors) 

Shaping  Filters  in  State  Equation 

x = Ax  + Bu  + Ew 

and  Noise  Power  W = cov(w) 

Workload  Expectations 

f 

c 

Attention  Allocation 

min  J 

f 'rfc.  = fc'f=,20 

C;  1 1 

SECTION  IV 
MODEL  DEVELOPMENT 


This  section  of  the  report  examines  some  of  the  experimental  data  reported  in 
the  literature  to  determine  appropriate  values  of  model  parameters  for  incorporating 
motion  cues  and  Visual  Meteorological  Condition  (VMC)  cues.  The  noise/signal  ratio 
for  the  vestibular  measurements  is  determined  from  roll  tracking  data  recorded  in 
Reference  17.  These  parameters  are  confirmed  by  applying  the  model  to  a VTOL  hover- 
ing task  described  in  Reference  23,  which  has  data  for  fixed-base,  angular  motion  only, 
and  angular  and  linear  motion  cues.  The  effect  of  VMC  cues  on  the  pilot  model  is 
determined  from  a comparison  of  the  pilot  describing  function  data  reported  in  Reference 
58.  The  predictions  of  the  model  with  combined  motion  and  VMC  cues  is  examined. 

4.1  MODEL  PARAMETERS  FOR  MOTION  CUES 

This  subsection  describes  the  inclusion  of  vestibular  measurements  in  the 
optimal  control  model  of  the  pilot  to  account  for  the  effect  of  motion  cues.  The  data 
base  compiled  by  Shirley  (Reference  17)  was  used  for  this  purpose. 


4.1.1 


DESCRIPTION  OF  THE  EXPERIMENT 


A block  diagram  of  the  system  used  by  Shirley  is  displayed  in  Figure  4,  which 
shows  the  various  parameters  of  the  simulation.  The  two  input  paths  to  the  human  are 
the  visual  system,  in  which  the  error  is  converted  by  a scale  factor,  K^,  to  distance 
on  a display,  and  the  motion  system,  in  which  the  error  is  converted  to  a roll  angle 
through  the  scale  factor  K j..  The  operator's  output  (in  cm)  is  converted  to  a voltage 
by  a gain,  K This  voltage  drives  the  plant  which  is  then  compared  with  the  input 

signal  (a  sum  of  10  sine  waves). 
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HYBRID 

COMPUTER 


ROLL  ANGLE  PROPORTIONAL  TO  SYSTEM  ERROR 

Kmot  (deg/volt)  CHOSEN  FOR  EACH  SET  OF  DYNAMICS 

Kdis/Kmot  = CONSTANT  THROUGHOUT  ENTIRE  EXPERIMENT 


Figure  4.  Block  Diagram  of  Shirley's  Roll  Tracking  Experiment. 

The  motion  system  gain  was  chosen  individually  for  each  set  of  dynamics  and 
resulted  in  reported  excursions  of  approximately  + 5°  for  easy-to-control  dynamics,  and 
+ 15°  for  difficult-to-control  dynamics.  The  ratio  of  display  gain  to  motion  gain  is 
kept  constant  throughout  the  experiment.  The  exact  value  of  the  motion  gain  for  each 
set  of  dynamics  is  not  available  from  the  data.  However,  an  estimate  of  this  gain  may 
be  obtained  from  the  RMS  tracking  error.  Table  2 presents  estimates  of  the  motion 
system  gain  for  both  easy  and  hard  dynamics  under  the  assumption  that  the  motion  limit 
(5°  or  15°)  is  the  2a  value.  Since  the  estimated  motion  system  gain  does  not  vary  sig- 
nificantly, we  have  assumed  a value  of  1 .6  deg/ volt  throughout.  Note  that  if  there 
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Table  2.  Calculation  of  Motion  System  Gain. 


.2  _ HI  volts  _se_c_  = ] ] >7  vo|ts2  for  ISI  = 1050 
89.6  sec 

rISE  = relative  Integral  Squared  Error 
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were  no  vestibular  or  visual  thresholds,  the  values  of  the  motion  and  the  display  gains 
would  not  be  a factor  in  the  optimal  control  model  of  the  human  operator. 


The  dynamics  of  the  motion  system  used  in  Shirley's  experiment  are  approxi- 
mated by  a 0.1  second  pure  time  delay  up  to  a bandwidth  of  10  radians  per  second. 
Following  Dillow  et  al.  (Reference  36),  this  time  delay  will  be  lumped  with  that  of 
the  human  operator,  a reasonable  assumption  when  one  views  the  location  of  the  time 
delay  in  the  block  diagram.  Thus  the  time  delay  parameter  of  the  human  operator  is 
set  at  0.3  second  instead  of  the  nominal  0.2  second. 

The  input  signal  for  Shirley's  experiments  was  a sum  of  ten  sine  waves  with  a 
10  to  1 shelf  at  2.5  radians  per  second  for  most  of  the  runs.  Dillow  et  al.  have  pro- 
posed approximating  this  by  a rational  spectrum  derived  from  a four-pole  filter  driven 
by  white  noise  with  covariance  W.  Dillow  et  al.  have  derived  a procedure  for  choosing 
an  optimum  break  frequency  for  a four-pole  filter  approximation;  for  these  conditions  it 
turns  out  to  be  2 rad/sec. 

To  find  the  value  of  W corresponding  to  the  known  value  of  mean  square 

2 

output  (in  volt  ),  we  note  that  the  spectrum  of  the  (approximate)  input  signal  is  given  by 
*..(s)  = — (21) 

(s  + 

where  idq  is  the  break  frequency  of  the  filter.  The  mean  square  value  of  the  filter 
output  is  given  by 

— iao 

$..(s)  ds  = residue  of  * . . (s)  at  s = -u>q  (??) 

ao 

To  find  the  residue  of  this  fourth  order  pole,  we  use  the  relation 
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which  provides 


.2  = 20W 

Co-)7  5 = "u,o 

Solving  this  expression  for  W in  terms  of  the  mean  square  output  provides 

W = — = 9601  volts^  sec,^  for  ISI  = 1050  (24) 

10  U 

Many  of  the  controlled  elements  considered  by  Shirley  were  second  order 
systems  (and  a pure  delay  of  0.1  second).  Thus  we  will  use  the  following  state  equations 
and  measurement  equations  for  the  roll  control  of  an  arbitrary  second-order  system  of  the 
form 

Y_(s)  = 1 (25) 

2 

S + Ps  + Of 

The  system  differential  equations  for  this  general  second-order  plant  are 


given  below: 

Shaping  Filter  (u)q  = 2.0) 


Controlled  Element 


x5  = X6 

x6  = -Q'X5-ex6  + YU 


Semicircular  Canal 

x7  = x8 

= -0.001852  x?  - 0.08888  xg  + Kmot'e 

= Kmotx3+Q,,<motx5+ei<motx6 

- 0.001852  x?  - 0.08888  xg  - «mot  yu 

where  e = Xj  - x^ 

Otolith 

The  specific  force  f in  g's  is 


f (in  g’s) 


Kmot 


0 (in  deg) 


where  i is  the  displacement  of  the  otolith  above  the  roll  axis.  The  state  equation 
becomes 


= -1 .289  • 10“5  x7  + 0.34785  xg  + 0.00696  • Kmot 


= 0.006960  Kmot  x3  + 0.006960  Kmof«  x5 

+ 0.00696  K tPx,-l  .289  • 10  5 x?  + 0.34785  xg 
mot  o 7 

- 0.006960  Kmot  Yu  (3°) 


Otolith  (g's) 

y = -1 .024  x0  + 2.024  f (in  g's) 

'oto  Y 

= t 2.024  K^,x,  - 2.024  K^,  ^3 

- 2.024  (l+^)»5 

- 2.024  K^0,  - 8 x6  - 1.024  x,  + 2.024  K|nof  - yu  (31) 

g 9 

hive  second-order  plants  were  chosen  from  Shirley’s  data  to  provide  a repre- 
sentative cross  section  of  controlled  elements.  These  are  shown  in  Table  3 with  the 
relative  Integral  Squared  Error  (HSE) and  standard  deviation  of  the  tracking  error  (ae).  Note 
"that  the  single  integrator  plant  has  been  approximated  by  a second-order  plant  with  high 
bandwidth  filter  to  avoid  the  analytical  (but  not  numerical)  problems  of  defining  the 
acceleration  input  to  the  vestibular  system.  The  system  matrices  for  the  differential 
equations  describing  the  controlled  element  and  measurements  are  shown  in  Tables  4a 
through  4e  for  these  five  plants.  Note  that  in  all  cases  the  E matrix  (noise  distribution 

matrix)  is  E ' = (0,0,0, 1 ,0,0, 0,0,0). 
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Table  3.  Roll  Control  Tasks  from  Shirley. 


Plant 

Vision 

rISE/ 

°e  (volts) 

Motion 

rISE/ 

°e  (volts) 

Vision/ 

Motion 

rISE/ 

°e  (volts) 

Shirley's 

Designation 

1 

.34 

.33 

.24 

Data  33 

s(s+4) 

1 .99 

1.96 

1 .67 

1 ~ 1 

.21 

.15 

.13 

Data  7-8 

S s (.02s  + 1) 

1 .56 

1 .32 

1 .36 

1 

1.09 

.57 

.44 

Data  19 

2 

s 

3.56 

2.58 

2.27 

15 

.39 

.32 

.30 

Data  47 

s2  + 3s  +15 

2.13 

1.94 

1 .87 

10 

.81 

.66 

.53 

Data  40 

s2  + 10 

3.08 

2.78 

2.49 

The  numerical  values  of  the  model  parameters  which  are  common  to  all  five 

controlled  elements  are.-  the  weighting  coefficients  on  displayed  quantities  (set  to  the 

inverse  square  of  the  threshold  values);  the  disturbance  covariance  (9601);  the  controller 

delay  (0.3  second,  which  includes  the  controlled  element  delay);  the  motion  noise/ 

2 

signal  ratio  (-25  dB);  and  the  vestibular  threshold  values  (0.19  deg/sec  , 0.005  g's). 

No  quadratic  weightings  were  placed  on  state  deviations  or  control  displacement  devia- 
tions, and  the  control-rate  weightings  were  adjusted  to  provide  a neuromuscular  time 
constant  of  0.1  second  (see  Table  5). 
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Table  4a.  System  Matrices  for  Shirley's  Plant  i/s  (s  + 4). 


0.003514 


Table  4b.  System  Matrices  for  Shirley's  Plant  1/s  (,02s  +•  1). 


0.1757 


Table  4d.  System  Matrices  for  Shirley's  Plant  1 5/(s  + 3s  + 15). 


-0.1670 


Table  4e.  System  Matrices  for  Shirley's  Plant  10/(s^  + 10). 


0.03514 


Table  5.  Control -Rate  Weighting  Coefficients  for  Shirley's  Cases. 

<tn  = 0,1  sec) 


Controlled 

Element 

Qr 

(sec^/in^) 

u 

max 

(in/sec) 

1 

0.012 

9.13 

s(s  + 4) 

1 

2.20 

0.67 

s (0.02s  + 1) 

1 

0.035 

5.35 

2 

s 

15 

2.476 

0.64 

s2  + 3s  + 15 

10 

2.468 

0.64 

s2  + 10 

The  only  remaining  unspecified  parameter  in  the  model  is  the  vestibular 

noise/signal  ratio,  which  is  the  one  under  investigation.  The  RMS  error  predictions  of 

the  model  are  compared  with  the  experimentally  observed  RMS  tracking  errors  in  Figure  5. 

It  can  be  seen  from  the  figure  that  three  of  the  five  controlled  elements  chosen  for 

analysis  give  a consistent  noise/signal  ratio  of  approximately  -18  dB  for  agreement 

between  the  model  and  the  experimental  data.  The  higher  predictions  of  the  l/s  plant, 

2 

and  the  lower  predictions  of  the  10/ (s  + 10)  plant  could  arise  from  the  effect  of 
different  levels  of  motor  noise  induced  by  motion  effects.  Since  the  l/s  plant  is 
generally  a very  easy  control  task  that  would  result  in  smooth  motion,  a lower  value 
of  signal  to  noise  ratio  would  be  expected  as  compared  to  the  highly  oscillatory  and 
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Figure  5.  RMS  Tracking  Error:  Comparison  of  Model  Predictions 
With  Shirley's  Experimental  Data. 


2 

more  abrupt  motions  resulting  in  the  control  of  lO/(s  + 10).  This  explanation  would  be 
appropriate  for  the  oscillatory  plant . but  is  not  entirely  satisfactory  for  the  integrator, 
since  the  -25  dB  value  is  appropriate  for  fixed-base  tracking  conditions. 

To  further  explore  the  predictions  of  the  model.  Figures  6a  through  6e 
depict  the  describing  functions  predicted  by  the  model  and  the  corresponding  experi- 
mental data  recorded  by  Shirley.  For  all  of  the  plants  except  1/s,  there  is  good  to 
excellent  agreement  between  the  experimental  results  and  the  analytical  predictions  in 
the  mid  to  high  frequency  range.  All  describing  functions  show  a characteristic  decrease 
in  amplitude  and  increase  in  phase  at  frequencies  below  approximately  0.3  rad/sec, 
which  is  consistent  with  a low  frequency  differentiator . On  the  other  hand,  the 
describing  function  for  l/s  has  the  same  basic  form,  but  is  obviously  operating  at  too 
low  a gain,  which  also  explains  the  model's  prediction  of  a larger  RMS  tracking  error 
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than  actually  observed.  The  low  frequency  discrepancies  between  the  experimental  data 
and  the  model  can  be  summarized  by  the  forward  loop  transfer  function  behavior;  the  ob- 
served forward  loop  transfer  function  would  appear  to  be  more  like  an  integrator  at  low 
frequencies,  whereas  the  model,  in  its  current  configuration,  appears  more  like  a high 
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cu  (rad/sec) 


Comparison  of  Model  Describing  Functions  Wi th  Shirl 
Experimental  Data.  (Y  = l/s  (s  + 4)  ) 


PHASE  (deg)  MAGNITUDE  (dB) 


0.1  1.0  10 

cu(rad/sec) 


Figure  6c.  Comparison  of  Model  Describing  Functions  With  Shirley's 
Experimental  Data.  (Y  = 1/s*) 


Gj(rad/sec| 


Figure  6e.  Comparison  of  Model  Describing  Functions  Wiih  Shirley's 
Experimental  Data.  (Y  = 10/(s^  + 10)) 


gain  at  low  frequencies.  This  can  be  shown  to  be  the  optimal  control  strategy  for  a 
first  order  correlated  input.  Thus,  the  low  frequency  behavior  of  the  model  appears 
to  be  an  inherent  characteristic  of  the  optimization  approach  and  the  modeling  of  the 
motor  noise  at  the  output  of  the  model. 

To  explore  the  appropriateness  of  the  cost  function,  the  optimal  control 
model  was  exercised  with  different  relative  weights  on  semicircular  canal  and  otolith 
measurements.  Our  hypothesis  was  that  an  increased  weight  on  the  otolith  measure- 
ment would  raise  the  pilot  gain  at  low  frequencies  because  the  low  frequency  otolith 
measurement  is  essentially  one  of  angle  with  respect  to  the  vertical.  However, 
increasing  the  weighting  on  the  otolith  measurement  in  the  cost  functional  had  neglig- 
ible effect  on  the  determination  of  the  feedback  gains,  hence,  the  control  strategy 
was  essentially  the  same  as  when  the  otolith  and  semicircular  canal  were  weighted 
equally.  Increasing  the  relative  importance  of  the  semicircular  canal  measurement 
provided  only  a more  pronounced  deviation  from  the  experimental  data  at  the  lower 
frequencies.  Consequently,  with  the  optimal  control  structure  of  the  pilot  model,  the 
low  frequency  open-loop  transfer  function  is  likely  to  appear  more  as  a gain  rather 
than  integrator,  and  this  may  differ  from  experimentally  observed  describing  functions 
in  some  instances. 

While  it  must  be  recognized  that  there  is  a deviation  in  the  describing 
function  predictions  at  low  frequency,  the  measurement  of  describing  functions  is  not 
reliable  at  these  frequencies,  and  nearly  the  same  control  performance  may  be  obtained 
in  this  portion  of  the  spectrum  if  the  human  operator's  gain  is  high  enough.  Thus,  it 
was  decided  to  continue  to  use  the  noise/signal  ratio  of  -18  dB  to  represent  the 
vestibular  system  remnant.  This  is  a good  prediction  of  RMS  performance  for  three 
plants,  and,  as  will  be  seem  subsequently,  it  gives  good  predictions  of  performance  in 
other  situations  as  well. 
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4.2  APPLICATION  OF  THE  MOTION  CUE  MODEL  TO  VTOL  HOVERING 

This  section  describes  the  application  of  the  optimal  control  model  of  the 


pilot  to  the  VTOL  hovering  task  described  in  Ringland,  Stapleford,  and  Magdaleno 
(Reference  23).  The  model  has  been  applied  to  the  following  motion  conditions: 

• Fixed-Base  (FB) 

• Moving-Base,  Angular  Motion  Only  (MBA) 

• Moving-Base,  Linear  and  Angular  Cues  (MBL) 

Three  tracking  tasks  were  conducted  simultaneously;  these  were: 

• Longitudinal 

• Lateral 

• Vertical 

For  the  longitudinal  and  lateral  tracking  tasks,  experimental  data  are  available  for  both 
"good"  and  "bad"  dynamics,  corresponding  to  the  degree  of  aerodynamic  damping;  a 
total  of  14  tracking  tasks  were  possible  in  these  conditions.  In  addition,  these  data 
provided  an  opportunity  to  compare  visual  scanning  behavior  with  and  without  motion 
cues  to  the  attention  allocation  among  instruments  as  predicted  by  the  model. 

4.2.1  DESCRIPTION  OF  THE  EXPERIMENT 

This  section  describes  the  experiment  used  in  the  VTOL  hovering  task  as 
described  in  Reference  23.  The  equations  of  motion  of  the  vehicle  are  given  by 

• Longitudinal 


s(s  - X )x  + g9  = -X  u 
' u'  a u g 

-M  sx  + s(s  - M )9  = M,  6 -Mu 
u ' q 6 e u g 


(32) 
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• Lateral 

s(s  - Yy)y  -gcp  = -Yyvg  J 

(33) 

-Lvsy  + s(s  - L^cp  = L^6a-Lvvg  J 

• Vertical 

s(s  - Z - T )z  = Z*  t>r  (34) 

' w w'  5 c 

c 

The  heading  of  the  simulator  was  held  constant  through  all  experimental  runs. 
The  disturbance  input  to  these  equations  were  generated  from  first-order  shaping  filters 
wi  th  a time  constant  of  1 second  (these  were  also  included  in  the  state  equations).  The 
state  equations  were  further  augmented  by  the  dynamics  of  the  simulator  drive  system, 
which  were  significant  in  the  longitudinal  direction.  As  shown  in  Reference  23,  the 
motion  system  can  be  approximated  by  a first-order  lag  with  a time  constant  of  0.1 
second.  The  aerodynamic  parameters  used  for  examining  the  various  motion  conditions 
are  shown  in  Table  6.  The  control  sensitivities  used  in  the  experimental  investigation 
are  given  in  Table  7;  we  chose  those  values  of  gain  corresponding  to  ’'good"  configu- 
rations for  each  of  three  tracking  tasks.  The  configurations  resulting  from  the  possible 
combinations  of  the  dynamics  are  shown  in  Table  8.  The  majority  of  experimental 
results  were  given  for  configurations  number  1 , 3,  4,  and  6 which  correspond  to  the 
four  possibilities  of  "good"  and  "bad"  dynamics  in  each  of  the  longitudinal  and  lateral 
axes.  These  are  the  systems  examined  with  the  model.  The  displays  used  by  the  pilot 
are  shown  schematically  in  Figure  7. 


Table  6.  Vehicle  Aerodynami  Parameters  (Reference  23) 


a) 

Longitudinal  Dynamics 

Fixed  Parameters: 

gM^  = 0.2  sec  ^ , 

X =-0.1 
u 

-1 

sec 

Variable  Parameters: 

M (sec  ') 

q 

Descriptor 

-4.0 

"Good" 

1.0 

"Mediocre" 

0 

"Bad" 

b) 

Lateral  Dynamics 

Fixed  Parameters: 

gL  = -0.2  sec  ^ , 

-< 

< 

ii 

i 

o 

-1 

sec 

Variable  Parameters: 

Lp(sec 

Descriptor 

-4.0 

"Good" 

-0.5 

"Bad" 

c) 

Vertical  Dynamics 

Variable  Parameters: 

Zw(sec  ')  Z^(sec” 

Descriptor 

-1.00  _3.00 

"Good" 

-1.00  -0 

"Mediocre" 

-0.25  0 

"Bad" 

Configuration 

Descriptor 

. 

Optimum 

Gain 

Subject 

RG 

EF 

GB 

"Good" 

0.066 

0.070 

0.043 

"Mediocre" 

0.047 

0.037 

0.037 

"Bad" 

0.042 

0.031 

0.031 

b)  Lateral  Task,  L ^ tad/sec 
a inch 


Configuration 

Descriptor 

Optimum 

Gain 

Subiect 

RG 

EF 

GB 

"Good" 

0.074 

0.089 

0.095 

0.053 

"Bad" 

0.040 



0.039 

0.095 

0.053 

c)  Vertical  Task,  Z;  Se— 
c inch 


Configuration 

Descriptor 

Optimum 

Gain 

Subject 

RG 

EF 

GB 

"Good" 

-4.3 

-4.91 

-4.91 

-4.91 

"Mediocre" 

-2.9 

-3.10 

-3.10 

-3.10 

"Bad" 

-2.2 

-0.77 

- 

-0.77 

-0.77 
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Table  8.  Vehicle  Configurations  (Reference  23) 


Configuration  No. 

Longitudinal 

Lateral 

Vertical 

* 

MMH 

r 

1 

"Good'' 

"Good" 

"Good" 

Oft! 

2 

"Mediocre" 

"Good" 

"Good" 

0 ft 

H 

"Good" 

"Good" 

11 

1 ■■{ 

■■■■■■■■■■ 

■■■■■■■■■ 

fil 

H 

"Good" 

"Bad" 

"Good" 

B 

5 

"Mediocre" 

"Bad" 

"Good" 

0 ft 

I * 

"Bad" 

"Bad" 

"Good" 

o ft ! 

i 

7 

"Good" 

"Good" 

"Mediocre 11 

0 ft 

8 

"Good" 

"Good" 

"Bad" 

0 ft 

9 

"Good" 

"Good" 

"Good" 

20  ft 

10 

"Mediocre" 

"Good" 

"Good" 

20  ft 

11 

■ — . 

"Bad" 

"Good" 

"Good" 

20  ft 

4.2.2  MODEL  PARAMETERS 

The  controlled  element  dynamics  for  the  control  task  are  decoupled  — thus 
the  complete  task  may  be  considered  as  the  simultaneous  control  of  three  independent 
systems.  The  control led-state  equations  were  derived  with  the  assignment  of  variables 

shown  in  Table  9.  The  resulting  system  matrices  are  shown  in  Tables  10  through  12. 

The  output  matrices  (C,  D)  were  derived  from  thte  visual  display  gains  and  the  models 

of  the  vestibular  sensors. 

The  display  gains,  weighting  coefficients,  and  indifference  thresholds  are 
summarized  in  Table  13.  The  weighting  coefficients  for  the  position  of  a display 
element  were  chosen  as  the  inverse  square  of  what  a pilot  would  tolerate  as  a maximum 
value.  The  indifference  threshold  was  chosen  in  the  same  manner,  generally  20  - 25 
percent  of  the  maximum  allowable  value.  Again,  it  is  relatively  straightforward  to 


Table  9.  State-Variable  Assignment  for  VTOL  Hovering  Task. 


State  Variable 

Longitudinal 

Lateral 

Vertical 

1 

Disturbance  Filter 

Disturbance  Filter 

Disturbance  Filter 

2 

X 

y 

z 

3 

X 

y 

z 

4 

Simulator  Lag* 

0 

Otolith 

(x) 

5 

0 

0 

6 

q 

see, 

7 

see, 

SCC2 

8 

SCC2 

Otol  ith 

9 

Otolith 

Included  only  for  MBL. 


64  - 


Table  10a.  RS&M  System  Matrices  - Longitudinal,  Good,  FB. 


Table  10c.  RS&M  System  Me  trices  - Longitudinal,  Good,  MBL. 
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Table  lOd.  RS&M  System  Matrices  - Longitudinal,  Bad,  FB 


Table  lOe.  RS&M  System  Matrices  - Longitudinal,  Bad,  MBA. 


-0.002479  0 0.002479  0 0 -0.00001289  0.3475 


Table  11a.  RS&M  System  Matrices  - Lateral,  Good,  FB. 


286.5 
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Table  lie.  RS&M  System  Matrices  - Lateral,  Good,  MBL. 


Table  lid.  RS&M  System  Matrices  - Lateral,  Bad,  FB. 


Table  lie.  RS&M  System  Matrices  - Lateral,  Bad,  MBA. 


0.00247  0 -0.002479  0 -0.1994  -0.00001289  0.3475 


Table  Ilf.  RS&M  System  Matrices  - Lateral,  Bad,  MBL. 


Table  13.  Display  Gains,  Weighting  Coefficients  ana  Indifference  Thresholds. 


Display 

Variables 

Display  Gain 

Weighting 

Coefficient 

1 ndifference 
Threshold 

1 in/9.1  ft 

1/(1  in)2 

.2  in 

u,v 

(1  in/sec)/(9.1  fps) 

(5  sec/1  in)2 

.2  in/5  sec 

e,<t> 

5 ( 1 80/r ) deg/ rad 

1/(10  deg)2 

2 deg 

o,<fi 

5 ( 180/tt)  de9/sec 
rad/sec 

(1  sec/10  deg)2 

2 deg/1  sec 

Z 

l .0 

1/(10  ft)2 

2 ft 

W 

l .0 

(5  sec/ 10  ft)2 

2 ft/5  sec 

Semicircular  ^ 
Canal 

1 .0 

0 

. 57  deg/sec2 

Otol  ith 

1.0 

0 

• 015  g 

obtain  weighting  coefficients  and  indifference  thresholds  on  position  of  a display  element 

(such  as  attitude),  but  it  is  more  difficult  to  assign  the  same  quantities  for  the  rates. 

In  Table  13  the  technique  suggested  in  Section  3 has  been  used  to  find  the  ratio  y / 

'max' 

^max  ' ',e*'  equivalent  time  constant  of  the  desired  response.  For  the  transla- 
tional variables  an  equivalent  time  constant  of  5 seconds  was  chosen,  whereas  the  same 
ratio  was  taken  as  1 second  for  the  attitude  control  task. 

Incorporating  both  visual  and  motion  cues  into  the  Kalman  filter  portion  of 
the  model  is  straightforward,  but  the  question  of  weighting  vestibular  "outputs"  must 
be  addressed.  In  the  motion-only  conditions  run  by  Shirley,  this  was  not  a problem  be- 
cause there  are  no  visual  terms  to  compete  with  the  vestibular  terms  in  a cost  functional. 
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One  attempt  was  made  at  choosing  weights  lor  the  vestibular  outputs,  in  a manner  simi- 

0 

lar  to  choosing  weights  for  display  element  rates.  If  a representative  ratio  of  0max/ 

• • • • • 

0 could  be  found,  then,  since  0 is  known,  0 could  be  calculated  and  the 
*max  rmax  max 

vestibular  weights  chosen  as  the  inverse  square  of  this  value.  This  would  allow  a 
"bootstrap"  determination  of  angle  rate  and  vestibular  weights  from  the  maximum  allow- 
able angle.  From  Dillow  et  al.  (Reference  36)  it  was  found  that  the  ratio  of  weights  on 

• • • 

roll  rate  and  acceleration  correspond  to  0 /0  =0.1  second.  This  relative 

' max  max 

weighting  was  tried  but  found  to  be  totally  unacceptable  because  of  too  much  weighting 
on  the  vestibular  outputs.  We  have  taken  the  opposite  tack  at  this  juncture  and  placed 
zero  weight  on  the  vestibular  measurements,  so  that  the  model  treats  the  vestibular 
signal  strictly  as  an  additional  measurement  used  in  the  control  feedback  loop,  but  not 
as  variables  to  be  minimized.  The  validity  of  this  approach  will  be  re-examined  later 
in  other  control  tasks. 

4.2.3  MODEL  RESULTS 

The  RMS  tracking  performance  for  each  of  the  three  tracking  tasks  was 
determined  under  the  various  motion  conditions  and  vehicle  dynamics  outlined  above. 
These  results  are  shown  in  Tables  14a  through  14c  for  different  levels  of  control 
workload,  f ' , the  fraction  of  attention  applied  to  the  control  task.  The  data  in  these 
tables  are  plotted  as  a function  of  workload  in  Figures  8 through  10.  The  characteris- 
tics to  note  in  these  plots  are  that  motion  cues  reduce  the  RMS  error  as  expected  in 
this  situation  of  veridical  motion,  and  the  “bad"  dynamics  generally  exhibit  more 
sensitivity  to  changes  in  workload  than  "good”  dynamics.  It  is  well  known  that  the 
slope  of  these  performance  versus  workload  curves  is  an  important  element  in  pilot 
opinion. 

At  each  level  of  control  attention,  fc,  the  allocation  between  the  position 
display  and  attitude  display  has  been  optimized.  The  ratio  of  the  control  attention 
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Table  14a.  Longitudinal  Performance  versus  Workload. 
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Table  14b.  Lateral  Performance  versus  Workload. 
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Table  14c.  Vertical  Performance  versus  Workload. 
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spent  on  the  position  and  attitude  displays  f ^/f  is  also  shown  in  Tables  14a  through 
14c,  which  indicate  that  bad  dynamics  generally  require  more  attention  to  the  attitude 
display  to  derive  the  necessary  information  for  stabilization,  whereas  with  good  dynam- 
ics, the  attention  can  be  devoted  to  the  ultimate  position  control  task. 

A prediction  of  overall  performance  for  the  three  motion  conditions  was 
obtained  as  follows.  Because  the  pilots  said  they  spent  very  little  time  on  the  altitude 
control  task,  we  used  the  vertical  performance  predictions  for  = 0 . 1 . To  obtain  the 
performance  predictions  for  the  lateral  and  longitudinal  axes,  the  fraction  of  attention 
on  lateral  and  longitudinal  axes  was  found  by  hand,  which  minimizes  the  total  RMS 
hovering  error,  sugject  to  the  constraint  that  the  attention  on  both  lateral  and  longitudi- 
nal axes  sum  to  0.9.  Results  of  our  previous  studies  have  indicated,  however,  that  the 
minimum  in  overall  RMS  position  is  generally  quite  flat  relative  to  the  allocation  of 
attention  between  longitudinal  and  lateral  tasks;  i.e.,  one  can  tradeoff  longitudinal 
versus  lateral  position  error  with  very  little  change  in  overall  RMS  position  error.  The 
only  requirement  for  this  insensitivity  is  that  enough  attention  be  applied  to  any  “bad” 
dynamics  to  be  beyond  the  knee  of  the  performance  versus  workload  curve.  The  results 
of  the  total  task  performance  predictions  are  shown  in  Figure  11,  which  compares  the 
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10.  Vertical  Performance  Versus  Workload 


Oh 


Figure  11.  Preliminary  Predictions  of  Total  Task  Performance  (Circles) 
Compared  to  Range  of  Means  of  Three  Subjects. 


RMS  values  of  important  state  variables  as  predicted  by  the  model  with  the  range  of 
mean  values  observed  for  three  subjects.  Thpre  is  excellent  agreement  between  the 
model  predictions  and  the  observed  range  of  subject  means,  both  on  a relative  basis 
(with  and  without  motion  cues)  and  on  an  absolute  basis. 

The  ratio  of  fractions  of  attention  to  the  attitude  and  position  displays  are 
shown  in  Table  15a.  To  determine  the  fraction  of  attention  allocated  to  the  position 
and  attitude  displays  for  each  of  the  configurations,  the  appropriate  entries  to  Table 
15a  were  added  to  provide  the  results  in  Table  15b,  which  shows  the  predicted  atten- 
tion allocation  as  a function  of  the  motion  and  configuration.  Note  that  for  any 
given  configuration,  there  is  relatively  little  change  in  attention  allocation  for  the 
various  motion  conditions,  whereas  there  is  a substantial  change  for  configurational 
differences. 

It  seems  logical  that  there  should  be  a correlation  between  the  fraction  of 
attention  predicted  by  the  extended  model,  and  eye  scanning  behavior.  Although 
several  attempts  have  been  made  to  link  the  optimal  control  model  results  with  overt 
eye  scanning  data,  only  qualitative  agreement  between  the  two  has  been  found. 

There  are  several  reasons  for  this:  pilots  are  trained  to  use  certain  scanning  behavior, 
and  it  may  be  difficult  to  modify  this  behavior  with  a change  in  control  task  difficulty; 
even  if  a pilot  is  viewing  an  instrument  foveally,  he  may  not  be  paying  attention  to 
it;  conversely,  if  a pilot  is  viewing  one  instrument  directly,  he  may  be  able  to  gather 
information  from  nearby  instruments. 

The  attention  allocations  of  Table  15b  are  correlated  with  the  measured  dwell 
fraction  averaged  over  all  motion  conditions,  as  quoted  in  Reference  23,  and  these 
results  are  shown  in  Figure  12.  For  the  dwell  fraction  observed  on  both  the  attitude 


Table  15a.  Fraction  of  Attention  Allocation  by 
Motion  and  Dynamics. 
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Table  15b.  Predicted  Attention  Allocation  by  Motion  and  Configuration. 


CONFIGURATION 


MOTION 

1 

3 

4 

6 

Average  for 
Motion 

fatt 

f 

pos 

fatt 

f 

pos 

fatt 

f 

pos 

fatt 

f 

pos 

fatt 

f 

pos 

FB 

.04 

.86 

.41 

.49 

.41 

.49 

.78 

ID 

.41 

.49 

MBL 

.04 

.86 

.41 

.49 

.41 

.49 

.78 

.12 

.41 

.49 

MBA 

.02 

.88 

.40 

.50 

.37 

.53 

.72 

.12 

.38 

.52 

Average  for 
Configuration 

.03 

.87 

.41 

.49 

.40 

.50 

.76 

.14 

display  and  the  position  display,  there  is  a fairly  linear  relationship  between  ths  dwell 
fraction  and  the  fraction  of  attention  derived  from  the  model.  The  eye  scanning  statis- 
tics are  relatively  insensitive  to  these  configurational  changes;  for  example,  the  dwell 
fraction  on  the  attitude  display  varies  between  0.5  and  0.6  over  the  range  of  configu- 
rations examined.  Thus  the  fraction  of  attention  derived  from  the  model  certainly  varies 
in  the  same  manner  as  the  observed  dwell  fraction  (and  the  relationship  appears  to  be 
fairly  linear),  although  there  is  a difference  between  the  attitude  and  position  displays.  A 
closer  correspondence  between  the  fraction  of  attention  and  dwell  fraction  might  be 
obtained  by  artificially  imposing  constraints  that  the  fraction  of  attention  on  the  attitude 
instrument  should  always  be  greater  than  (say)0.5toaccountforpilotscanning  habits.  The 
attention  allocation,  as  used  in  the  present  model,  is  a normative  description  (i.e., 
how  the  pilot  should  allocate  his  attention).  Thus,  the  examination  of  more  eye  move- 
ment data  may  indicate  that  pilot  attention  will  have  to  be  bounded  to  reflect  pilot 
limitations.  This  result  would  be  consistent  with  the  other  constraints  currently  in  use 
with  the  optimal  control  model:  the  pilot  minimizes  the  quadratic  cost  functional  sub- 
ject to  constraints  on  his  inherent  capabilities. 

4.3  MODEL  PARAMETERS  FOR  VMC  CUES 

The  characteristics  of  the  human  operator  when  utilizing  VMC  cues  have  been 
qualitatively  described  by  use  of  larger  attitude  angles  to  maintain  a tighter  control  on 
position  (Reference  23).  Unfortunately,  this  reference  contains  no  substantial  data  which 
could  be  used  to  explore  the  various  manifestations  of  VMC  cues  within  the  model.  How- 
ever, Reference  25  provides  a fairly  complete  set  of  VMC/IMC  experimental  differences 
in  which  the  effects  of  VMC  cues  can  be  summarized  as  allowing  the  pilot  to  operate  at 
a higher  gain  than  under  IMC  cues.  As  mentioned  in  Section  3,5,  this  behavior  can  be 
introduced  into  the  model  by  three  methods:  decreasing  the  noise/signal  ratio  on  a tti  — 
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tude  rate  (corresponding  to  increased  acuity  of  roll  rate);  decreasing  the  threshold  on 
roll  rate  (again  due  to  increased  acuity  of  roll  rate);  and  decreasing  the  quadratic 
weighting  penalty  for  roll  rate  (thus  allowing  the  system  to  operate  at  a higher  loop 
gain).  The  first  two  methods  are  more  physiologically  based,  whereas  the  last  technique 
(similarly  for  increased  roll  angles  in  the  VTOL  hovering  task)  are  based  on  differences 
in  pilot  behavior. 

4.3.1  EXPERIMENTAL  DESCRIPTION 

To  explore  the  differences  in  angular  rate  threshold  (and  effectively,  noise/ 
signal  ratio)  and  the  effect  of  attitude  rate  on  the  model,  we  chose  a fixed-base  track- 
ing task  for  which  pilot  describing  function  data  have  been  obtained  by  Newell  and 
Smith  (Reference  25).  In  both  cases  the  controlled-element  dynam:cs  are  K/s(s  - 1), 
with  the  primary  difference  being  in  the  form  of  the  display:  a conventional  attitude 
indicator  (IMC  cues)  as  shown  in  Figure  13;  and  a wide  field  contract  analog  display 
(VMC  condition)  as  shown  in  Figure  14.  The  command  input  for  each  case  was  a sum 
of  sine  waves  with  a rectangular  spectrum  (shelf  frequency  at  2.5  radians/second). 

The  system  matrices  for  both  the  IMC  and  VMC  tracking  tasks  are  shown  in  Table  16a. 

The  IMC  tracking  condition  is  the  conventional  use  of  the  optimal  model  of 
the  human  operator.  The  model  parameters  to  predict  the  results  of  this  IMC  tracking 
task  are  shown  in  Table  16b.  The  maximum  allowable  roll  deviation  was  based  on  the 
display  by  the  attitude  indicator,  and  the  indifference  threshold  was  chosen  at  2 
degrees  because  of  the  layout  of  the  display.  Tfe  weights  and  thresholds  for  the  roll 
rate  were  chosen  on  the  basis  of  a 1 second  time  constant  response  for  this  attitude 
control  task.  The  control  rate  weighting  was  chosen  to  obtain  a neuromuscular  time 
constant  of  approximately  0.1  second. 
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’igure  14.  Contact-Analog  Display  for  VMC  Experiments  (Reference  25). 


Table  l6a„  System  Matrices  for  Newell-Smith  K/s(s  - 1). 
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Table  16b.  Model  Parameters  for  IMC  Tracking. 


Input  Variable 

Weight 

Threshold 

0 (deg) 

1/(5  deg)2 

2 deg 

0 (deg) 

(1  sec/5  deg)2 

2 deg/sec 

To  evaluate  the  effect  on  the  model  predictions  of  changing  roll  rate  thresh- 
olds and  roll  rate  weightings,  the  describing  functions  were  calculated  for  thresholds  of 
2 deg/sec,  0.7  deg/sec  and  0.2  deg/sec,  and  for  maximum  roll  rates  of  5 deg/sec, 

10  deg/sec  and  50  deg/sec.  In  each  one  of  these  conditions  there  is  an  overlap  with 
the  nominal  IMC  tracking  case,  as  shown  in  Table  16b.  The  effect  on  the  describing 
function  of  changing  roll  rate  threshold  is  shown  in  Figure  ]5,  along  with  the  experi- 
mentally measured  data  for  the  IMC  tracking  task.  It  is  obvious  that,  under  these 
conditions,  changing  the  roll  rate  threshold  has  almost  no  effect  on  the  describing 
function,  primarily  because  the  RMS  roll  rate  is  substantially  above  the  threshold 
value  even  at  the  highest  threshold  (2  degrees  per  second).  Essentially  then,  all  of 
these  describing  functions  are  predictions  of  the  describing  function  for  the  IMC 
tracking  conditions,  and  it  is  seen  that  there  is  good  agreement  between  the  experiment 
and  the  predictions.  Note  that  the  experimental  data  as  well  as  the  model  show  a 
fall-off  in  gain  at  the  lower  frequencies  under  these  conditions. 

The  effect  on  the  predicted  describing  function  for  changing  maximum  roll 
rate  is  shown  in  Figure  16.  Here  we  see  a substantial  alteration  of  the  describing  func- 
tion with  the  change  of  this  parameters,  indicating  that  it  is  a more  appropriate 
representation  of  VMC  tracking  cues  under  these  conditions.  The  experimental  data  in 
this  figure  are  the  describing  functions  determined  from  the  wide-field  contact  analog 
display  shown  in  Figure  14.  The  comparison  of  the  experimental  data  in  Figures  15  and 
16  indicate  that  when  operating  with  the  contact  analog  display,  the  pilot  can  obtain 
substantially  higher  gains,  especially  at  lower  frequencies.  In  this  particular  case,  the 
maximum  roll  rate  of  50  deg/sec  in  the  quadratic  penalty  provides  an  excellent  descrip- 
tion of  the  experimental  data. 
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Figure  15.  Effect  of  Attitude  Rate  Figure  16. 

Threshold  on  Pilot 
Describing  Function 
(Reference  25). 


Effect  of  Attitude  Rate 
Cost  Weighting  on  Pilot 
Describing  Function 
(Reference  25). 


For  the  roll  control  tracking  task  described  here,  the  effect  of  a wide  view 
contact  analog  display  (an  approximation  to  VMC  conditions)  seems  to  be  well  repre- 
sented by  an  increase  in  the  maximum  roll  rate,  rather  than  a decrease  in  roll  rate 
thresholds.  This  is  true  under  these  conditions  because  the  RMS  roll  rates  are  substan- 
tially above  that  threshold  even  in  IMC,  and  reduction  of  the  rate  thresholds 
does  not  improve  tracking  ability.  The  increase  in  pilot  gain  under  VMC  may 
arise  from  a desire  for  smaller  roll  angles  or  higher  roll  rates;  there  are  two  ways  to 
differentiate  between  these  conditions  in  terms  of  model  parameters.  In  the  first  case, 
if  the  RMS  tracking  error  is  at  or  slightly  above  the  indifference  threshold,  then  a 
reduction  in  roll  angular  error  would  not  seem  to  be  an  appropriate  goal,  thus  an 
increase  in  maximum  roll  rate  is  suggested.  The  second  case  determines  the  goal  by 
the  scaling  of  the  cost  functional.  Decreasing  the  roll-rate  weighting  coefficient  has 
no  appreciable  effect  on  the  control-rate  weighting  (maximum  control  in  inches  per 
second),  whereas  increasing  the  roll  angle  weighting  coefficient  is  likely  to  increase  the 
control  rate  weighting  coefficient  by  a commensurate  amount,  thus  implying  an  unreas- 
onable maximum  control  rate  to  obtain  the  desired  neuromuscular  time  constant  of  0.1 
second. 

4.4  APPLICATION  OF  THE  MODEL  WITH  MOTION  AND  VMC  CUES 

The  describing  functions  predicted  by  the  model  were  compared  with  experi- 
mental data  recorded  in  Newell  and  Smith  (Reference  25)  for  a roll  tracking  task  with 
IMC  and  VMC  cues.  Motion  cues  are  present  in  both  cases,  but  the  motion  effects  enter 
in  different  ways,  as  described  below.  Our  approach  to  modeling  the  VMC  cues  in  a 
tracking  task  is  to  decrease  the  maximum  roll  rate  as  described  in  Section  4.3. 
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4.4.1  DESCRIPTION  OF  THE  EXPERIMENT 

The  block  diagrams  of  the  two  system  configurations  used  in  the  study  are 
shown  in  Figure  17.  Configuration  C corresponds  to  Newell  and  Smith's  configuration 
A-2*  and  configuration  D corresponds  to  their  configuration  A-2.  The  mnemonic  labels 
C and  D are  used  to  correspond  to  Command  input  and  Disturbance  input,  respectively. 
Because  of  the  preshaping  of  the  (sum  of  sine  waves)  input  in  configuration  D,  both 
configurations  are  identical  in  fixed-base  experiments.  During  the  T-33  flight  experi- 
ments, configuration  C was  used  with  IMC  cues,  and  configuration  D was  used  with 
VMC  cues.  With  configuration  C in  flight,  the  roll  angle  error  indicated  on  the  display 
does  not  correspond  to  the  actual  aircraft  bank  angle.  The  motion  cues  transmit  the  air- 
craft motion  but  not  the  input  command  portion  of  the  displayed  roll  error,  similar  to  the 
experiments  conducted  by  Junker  and  Replogle  (Reference  39). 

4.4.2  MODEL  PARAMETERS 

To  approximate  the  shelf  spectrum  of  the  sum-of-sine-waves  input,  we  have 
used  a fourth-order  Butterworth  filter  which  is  maximally  flat.  The  state  equations  for 
the  shaping  filter  are 


*1  *2 


*2  = x3 


x3  = x4 


*a  = -a). 


x,  - 2.613u)  x~  - 3.414a,  x^  - 2.613<u  x . + tu  w 

O Z O o O H o 


The  only  information  given  in  Newell  and  Smith  is  that  the  maximum  amplitude  of  the 
command  was  8°.  If  we  assume  this  to  be  the  3ct  limit  on  a Gaussian  distribution,  then 
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Figure  17.  Block  Diagram  of  Control  Tasks  (NASA  TN  D-500] 


the  variance  of  the  white  driving  noise  in  Equation  (35)  should  be  chosen  such  that 
the  RMS  command  level  (x^)  is  2.667°. 

The  fifth  and  sixth  state  variables  represent  the  equations  for  the  controlled 

elements 


i 

‘i 


(36) 


The  righthand  sides  of  these  equations  will  be  derived  for  both  configurations  below. 

Equations  for  the  vestibular  system  are  the  same  as  described  in  Section  3.4 
and  are  given  by 

x7  = x8 

xg  = -0.001852  Xj  - 0.08888  Xg  + od  (deg/sec  ) (37) 

i9  = 0.2  x^,  + 0.2f(g's) 


The  output  of  the  vestibular  systems,  i.e.,  the  measurements  are  given  by 


y = -1 .289  . 10"5  x,  + 0.34785  xQ  + 0.00696  i 
see  / o 

yofo  = -1 .024x9  + 2.024  f 


(38) 


For  configuration  C,  the  only  input  to  the  controlled  element  is  from  the 
pilot's  control  stick  deflection  (here  denoted  u).  The  differential  equation  for 
roll  rate  in  flight  becomes  (Reference  25) 

= w = -2.857  0^  + 45u  (39) 

The  rate  of  change  of  roll  rate  for  configuration  D is  more  complex  because  of  the 
shaped  input  to  the  controlled  element.  This  equation  is  given  by 
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4>a  - 


x6 


-2.857  x6+  5 


u + 2.857  i 


45 


(40) 


x2  ” x3  ” 2*857  x6  + 45u 


where  the  second  portion  of  the  bracketed  expression  is  the  output  of  the  shoping  filter. 

The  specific  force  is  an  important  element  in  the  otolith  equations  and  is 
derived  from  three  sources:  gravitational  forces,  aerodynamic  forces,  and  constraint 
forces  from  linear  acceleration  due  to  roll  accelerations  with  a displacement  from  the 
roll  axis.  In  general  the  lateral  component  of  specific  force  (in  g's)  is 

,2 


r 2"  /*  M V‘  . £0A 

f = (sin  0A  - — cos  0.  

360  gR  g 


(41) 


where  V is  the  linear  velocity  and  R is  the  turn  radius.  If  there  is  no  translation  of  the 
center  of  gravity  during  roll  (R  = oo,  lateral  component  of  lift  is  cancelled  by  side 
force)  then  for  small  roll  angles  the  specific  force  (in  g's)  is 


f = ill 

360 


(•*-«*) 


(42) 


where  l is  the  distance  from  the  roll  axis  to  the  otoliths.  If  the  turn  is  aerodynamically 
coordinated,  then  the  lateral  component  of  gravi tational  specific  force  is  counteracted 
by  centripetal  acceleration  in  the  turn,  thus 


r _ 2 TT  Si  \ •' 

f 0 

360  g 


(43) 


The  former  condition  will  be  assumed  in  further  development  of  the  equations,  with  the 
realization  that  the  second  condition  can  be  derived  by  eliminating  the  0^  (x^)  term 
below.  Thus  the  specific  force  for  configuration  C is  (for  l = 4 ft): 

f = 0.01 745  x5  + 0. 00620  x6  - 0.09765u  (44) 
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The  corresponding  expression  for  the  specific  force  in  configuration  D is  given  by 


f = 0.00620  x2  + 0.002170  x3  + 0.01745  x5 
+ 0.00620  x6  - 0.09765u 


(45) 


It  is  now  possible  to  see  the  major  differences  between  configurations  C and 
D during  moving  base  conditions.  Equations  (39)  and  (40)  describe  the  input  to 
the  semicircular  canals  during  the  motion  condition.  In  configuration  D,  since  the 
attitude  of  the  aircraft  is  displayed  directly,  the  roll  angle  is  the  tracking  error,  so 
that  the  semicircular  canals  are  receiving  the  second  derivative  of  the  error  signal.  In 
configuration  C the  error  signal  is  (i  - 0^),  so  that  the  vestibular  system  is  not  receiving 
the  second  derivative  of  the  error  signal.  However,  the  motion  felt  is  due  to  the  pilot's 
control  inputs,  hence,  it  might  be  of  value  in  separating  command  from  response  in  the 
visually  displayed  error. 

Tables  17  and  18  show  the  system  matrices  used  in  the  study.  Table 
19  shows  the  assumed  values  of  weights  and  thresholds  derived  from  the  description 
of  the  experimental  setup  in  Newell  and  Smith.  Note  that  the  only  difference  between 
IMC  and  VMC  conditions  is  the  decreased  weight  on  roll  rate  in  VMC  conditions.  The 
thresholds  seem  to  be  reasonable  in  VMC  conditions  with  a well  defined  horizon  as  a 
reference  for  the  grease  pencil  stripes  on  the  windscreen  as  described  by  Newell  and 
Smith.  In  reality,  we  do  not  know  how  well  defined  the  horizon  was  in  these  tests,  so 
these  indifference  thresholds  may  be  optimistic  values. 

4.4.3  MODEL  RESULTS 

Figures  18  and  19  show  the  predictions  of  the  extended  pilot  model  with 
motion  cues  and  VMC  cues  along  with  the  experimental  describing  functions  measured 
in  flight.  For  these  describing  functions,  the  control-rate  weighting  was  based  on 
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Table  17.  System  Matrices  for  Newell-Smith  Configuration  C. 


0.3132  -0.1976) 


r 


Table  19.  Display  Weights  and  Thresholds. 


Measurement 

Weight 

Threshold 

IMC 

VMC 

IMC 

VMC 

0 (deg) 

0.04 

0.04 

2 

2 

0 (deg/sec) 

0.04 

0.0004 

2 

2 

Semicircular  Canal 
(deg/sec2) 

0 

0 

0.57 

0.57 

Otolith  (g) 

0 

0 

0.015 

0.015 

maximum  physical  rate  rather  than  the  neuromuscular  time  constant  of  0.1  second. 

Other  than  the  previously  described  motion  differences  which  have  been  incorporated 
in  the  model,  the  primary  difference  between  the  IMC  and  VMC  model  parameters  is 
the  decreased  weight  (i.e.,  increased  maximum  value)  on  roll  rate  for  VMC.  As  can 
be  seen  from  these  figures,  the  experimental  data  agree  well  with  the  predictions  in  the 
mid  to  high  frequency  range.  The  discrepancy  in  the  low  frequency  range  is  analogous 
to  the  deviations  observed  in  Shirley's  roll  tracking  data,  i.e.,  the  model  forward  loop 
transfer  function  behaves  as  a gain  at  low  frequencies  instead  of  an  integrator  as  implied 
by  the  experimental  data. 

This  point  was  pursued  further  by  incorporating  weights  on  the  otolith  mea- 
surement, on  the  premise  that  an  increased  weight  on  otolith  deviation  would  be 
analogous  to  increasing  the  weight  on  roll  error  at  low  frequencies.  However,  the 
increased  otolith  weights  produced  no  significant  change,  and  only  minor  changes 
occurred  when  weights  on  the  semicircular  canal  measurements  were  incorporated. 

Thus  it  appears  that  this  low  frequency  behavior  of  the  describing  function  cannot  be 
resolved  by  simply  adjusting  visual/vestibular  weighting,  but  appears  to  be  a charac- 
teristic of  the  optimal  control  model  as  it  is  currently  formulated. 
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In  spite  of  this  problem,  the  agreement  between  the  model  and  the  analytical 


predictions  is  quite  good  over  the  frequency  range  of  most  interest,  and  as  has  been 
seen  for  the  VTOL  hovering  task,  the  model  is  quite  capable  of  predicting  the  RMS 
performance  and  scanning  behavior  associated  with  motion  cue  effects. 


SECTION  V 


I 


1 


MODEL  APPLICATION 

In  this  section,  the  model,  with  the  parameters  determined  from  the  develop- 
ment phase,  is  applied  to  two  problems  for  which  experimental  data  exist  with  motion 
cues.  The  first  experiment  was  conducted  by  Stapleford,  Peters  and  Alex  (Reference  15) 
and  consisted  of  a 2 degree-of-freedom  roll  control  task  in  the  NASA  6 degree-of-freedom 
VTOL  simulator.  The  second  set  of  experiments  consisted  of  data  recorded  during  a simu- 
lated CTOL  approach  to  landing,  in  which  both  VMC  and  IMC  cues  were  available,  but 
in  which  the  predominant  cues  were  obtained  from  the  instruments. 

5.1  2 DEGREE-OF-FREEDOM  ROLL  CONTROL  TASK 

5.1.1  EXPERIMENT  DESCRIPTION 

The  experimental  description  in  Reference  15  is  quite  complete,  and  will  be 
described  here  only  as  it  pertains  to  that  portion  of  the  experiment  on  which  the  model 
was  used.  Controlled  element  dynamics  were  always  of  the  form 


K 

c 

s(s  + a) 


(46) 


where  a = 0,  1,  or  10  sec  \ The  controller  was  a low  inertia,  spring  restrained,  side 
stick  with  ±15  degrees  of  travel . It  had  a breakout  torque  of  approximately  1 . 5 in-lbs 
and  a gradient  of  approximately  0.32  in-lb/deg.  The  pilot-selected  gains  (Kc)  for  each 
of  the  four  subjects  are  shown  in  Table  20. 


The  overall  block  diagram  of  the  simulation  is  shown  in  Figure  20.  The 
disturbance  input  (d)  was  composed  of  10  sine  waves  having  a rectangular  shelf  spectrum, 
with  a break  frequency  at  1 .44  rad/sec  and  an  rms  input  of  0.4  deg. 
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Table  20.  Pilot-jelected  Gains. 


: 


a 

(sec-1) 

Subject 

GB 

GC 

RG 

MJ 

0 

6.5 

5 

5 

7.5 

1 

2.5 

4 

5 

10 

10 

65 

15 

25 

75 

Gain,  Kc,  in  units  of  sec 


The  display  consisted  of  the  roll  angle  error  shown  to  the  pilot  on  a 5-inch 
attitude  indicator,  and  was  inertia  1 1 y stabilized. 

The  roll  angle  washout  circuit  has  the  transfer  function 


0.  s 

-J-  = Y = 

0>  7 s + P 


(47) 


0 


where  p^  is  an  experimental  parameter  having  values  0.5,  1,  and  2 sec  \ To  simulate 
the  low  side  force  characteristics  of  a VTOL  in  hover,  the  following  transfer  function 
between  lateral  position  and  roll  angle  was  used 


y.  + bs  + g 

^ = y = z 5 

y <*v 


(48) 


where  is  the  simulated  pilot  displacement  along  the  z axis,  = 1 .5  ft  is  the  actual 
displacement  in  the  simulator,  b is  an  experimental  parameter  (usually  0),  and  p^  is  an 
experimental  parameter  having  values  0.5,  1,  or  2 sec  ^ . 

The  simulator  dynamics  were  tested  and  found  to  be  closely  approximated  by 
a first-order  lag  in  the  roll  dynamics: 
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INPUT 


Schematic  Diagram  for  2-DOF  Roll  Control  Task 


(49) 


0 -TS 

s • e 


0.  Ts  + 1 

i 

where  0 = actual  roll  angle  of  simulator  cab 

s 

0.  = commanded  roll  angle 
i 

t = 0.04  sec 
T = 1/12  sec 

The  lateral  position  dynamics  show  run  to  run  variability  for  frequencies 
greater  than  3 radians  per  second  (it  was  suggested  that  the  acceleration  elements  were 
being  reached).  An  approximation  to  the  position  dynamics  of  the  simulator  drive  system 
is  given  by 

2 


s 


yi  s2  + 2csu>ss  + u)2 


(50) 


where  y$  = actual  lateral  position  of  simulator  cab 

y.  = commanded  lateral  position 

C = 0.7 
ws 

(D  = 5 rad/sec 
s 

The  state  variable  assignment  for  the  optimal  control  pilot  model  is  shown  in 
Table  21,  where  again  the  rectangular  disturbance  spectrum  has  been  modeled  by  a 
fourth-order  Butterworth  filter.  The  displayed  variables,  their  weighting  coefficients  in 
the  quadratic  cost  functional , and  the  threshold  values  are  shown  in  Table  22.  The 
system  matrices  for  the  optimal  control  model  are  shown  in  Table  23  For  the  general 


case. 
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Table  21.  Stale  Variable  Assignment  for  2 Degree- 
of-Freedom  Roll  Control  Task. 


State 

Variable 

Simulation  Component 

*1  ' x4 

disturbance  shaping  filter 

x5'  x6 

controlled  element 

x7'  x8 

roll  ongie  dynamics  0/$. 

o 

X 

cs 

X 

Y\/<t> 

CM 

>r 

V 

X 

y/y\ 

y13'  x14 

semicircular  canal  dynamics 

X1 5 

. 

otolith  dynamics 

Table  22.  Displayed  Variables. 


Display 

Variable 

Weighting 

Coefficient 

Threshold 

0 

1/(5  deg)2 

2 deg 

0 

(1  sec/5  deg)2 

2 deg/sec 

see 

0 

0.57  deg/sec2 

OTO 

0 

0.01 5 g 

5.1.2  MODEL  RESULTS 


Figures  21  through  23  compare  the  describing  functions  for  moving-base 
conditions  and  fixed-base  conditions  obtained  with  the  optimal  control  model.  Also 
shown  are  the  exper  nental  results  obtained  by  Stapleford,  Peters  and  Alex  with  sub- 
ject RG  in  the  nor  ncil  condition  (P^=  P^  = .5  sec”',  ^ = b = 0).  These  figures  indi- 
cate that  both  the  model  and  the  experimental  data  show  only  small  differences  be- 


tween the  fixed-base  and  moving-base  conditions.  The  model  and  the  experimental 
data  both  agree  that  the  closed  loop  describing  function  gain  is  slightly  higher  for 
moving-base  conditions.  Above  approximately  1 rad/sec  the  model  agrees  very  well 
with  the  experimental  data;  however,  as  discussed  in  Section  5,  the  model  predicts 
somewhat  lower  gain  and  less  phase  lag  at  the  lower  frequencies. 


Figure  24  compares  the  mean-squared  errors  for  fixed-base  and  moving-base 
conditions  with  the  experimental  results  for  three  controlled  plants.  The  errors  have 
been  normalized  by  the  mean-squared  disturbance  input.  The  model  predictions  show 
the  same  trend  as  the  experimental  results.  However,  the  model  predicts  consistently 
large  errors  with  motion  and  somewhat  smaller  errors  for  the  fixed-base  condition. 


To  show  the  effects  of  changing  the  disturbance  input  spectrum,  the  Butter- 
worth  filter  was  replaced  by  a fourth-order  Dillow  filter  for  the  plant  K c/s(s+10)  with 
the  rms  input  level  unchanged.  The  closed  loop  describing  functions  for  fixed-base 
end  moving-base  are  compared  with  the  experimental  data  in  Figure  25.  Note  that 
changing  the  input  filter  has  as  much  effect  as  the  difference  between  moving-base  and 
fixed-base  conditions. 


The  effects  of  varying  the  simulator  dynamics  and  thereby  the  vestibular 

2 

inputs  were  examined  for  controlled  plant  K . Figure  26  shows  the  closed  loop 
describing  function  for  variations  in  the  lateral  position  washout  parameter  p^.  This 
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Figure  21 . 


Comparison  of  Fixed-Base  and  Moving-Base  Results  for  K/s(s+l) 
(Subject  RG). 
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parameter  changes  just  the  linear  acceleralion  cues  for  the  pilot;  varying  from  0.5 

sec-'  to  sec-'  showed  no  significant  effects  on  the  describing  function  nor  the  RMS 
results.  This  conclusion  was  confirmed  by  the  experimental  data  in  Reference  15. 

Similarly,  varying  the  roll  washout  parameter,  p^,  from  0.5  sec  to  2 sec  had 
very  little  effect  on  the  predicted  results,  as  shown  in  Figure  27.  This  variation 
alters  both  the  rotary  motion  cues  and  the  linear  acceleration  cues.  Again  the  model 
predictions  agree  with  the  observed  experimental  results. 

5.2  PILOT  PERFORMANCE  IN  LOW  VISIBILITY  APPROACH 

In  this  section  the  extended  pilot  model  is  used  to  predict  system  performance 
during  a low  visibility  approach.  The  experimental  data  was  generated  during  two 
phases  of  a low  visibility  approach  simulation  reported  in  detail  in  Reference  58 
(Phase  1)  and  Reference  59  (Phase  2). 

5.2.1  EXPERIMENT  DESCRIPTION 

The  aircraft  for  this  simulation  was  the  military  C-135B  transport  which  is 
quite  similar  to  the  commercial  Boeing  707  (Reference  60).  The  configuration  for  the 
experiment  was  fixed  at  a weight  of  160,000  lbs,  flaps  setting  at  50°,  and  gear  down. 
The  experimental  run  began  from  a trimmed  flight  condition  at  956  feet  altitude  on  a 3° 
glide  slope.  The  aerodynamic  coefficients  were  updated  as  a function  of  the  current 
flight  condition,  but  did  not  change  significantly. 

A conventional  flight  control  system  was  used  in  the  simulation:  the  controls 
included  wheel  (aileron),  column  (elevator),  and  rudders  (yaw)  and  four  throttle  settings. 
Figure  28  shows  the  instrument  panel  layout  used  in  the  simulation. 

The  cockpit  motion  system  had  3 degrees  of  freedom  — pitch,  roll,  and  heave 
within  the  following  limits: 
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Figure  27.  Effect  of  Roll  Washout  Variations. 
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vertical  position:  +13  inches 
pitch  angle:  +14°,  -6° 
roll  angle:  +9° 
vertical  acceleration:  + .8g 

2 

initial  angular  acceleration  in  roll:  50  deg/sec 

2 

initial  angular  acceleration  in  pitch:  50  deg/sec 
5.2.2  MODEL  DESCRIPTION 

The  model  was  applied  to  the  longitudinal  control  task  of  the  C-135B 
approach.  The  statistical  results  in  Reference  59  divide  the  approach  into  two  seg- 
ments based  on  the  range  from  touchdown: 

1 . Approach:  10,000  ft  - 5,000  ft 

2.  Final  approach:  5,000  ft-  1,500  ft 

Because  the  ASI  optimal  control  pilot  model  is  a stationary  one,  a frozen  point  analysis 
was  performed  for  the  midway  point  in  the  initial  approach  segment;  i.e.,  7,500  ft  from 
touchdown . 

The  model  disturbance  input  was  chosen  to  represent  the  Mil  Spec  8785B 
random  wind  described  in  Reference  61 . The  rms  gust  disturbances  were  2.5  fps, 

1 .65  fps,  and  ,33°/sec  for  the  longitudinal,  vertical  and  pitch  inputs  respectively 
(Reference  59).  The  dynamic  equations  of  motion  were  written  in  the  aircraft  body 
axis  reference  frame,  assuming  stabilized  flight  on  a 3°  glide  slope.  At  7500  ft  from 
touchdown,  this  corresponds  to  an  altitude  of  393  ft;  however,  since  the  stability 
derivatives  do  not  change  significantly,  the  tabulated  values  from  Reference  58  were 
used  for  an  altituae  of  500  ft.  These  data  are  summarized  in  Table  24.  The  simulator 
dynamics  were  approximated  by  first-order  lags  of  0.3  sec  in  heave  and  0.5  sec  in 
pitch  (Reference  60).  The  state  variable  assignments  for  the  model  are  summarized  in 
Table  25,  and  the  linear  equations  of  motion  are  shown  in  Table  26. 
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VT  = 233  ft/sec 


u = 233  ft/sec 
o 

w = 6.358  ft/sec 
o ' 

Yq  = "3  de9 

0 = -1 .439  deg 

o a 

h = 393  ft 

o 

R = 7500  ft 
o 

STABILITY  DERIVATIVES 


X 

u 

= -0.013632  sec-1 

Zu 

= -0.2553  sec  1 

M 

u 

= 0.000  3245  rad/sec/ft 

X 

w 

= 0.09922  sec"1 

Zw 

= -0.6274  sec"1 

M 

w 

= 0.004662  rod/sec/ft 

x. 

w 

= 0.0 

z- 

w 

= 0.0 

M- 

w 

= -0.001216  rad/ft 

X 

q 

= 0.0  ft/sec 

Z 

q 

= 0.0  ft/sec 

M 

q 

= -0.7167  sec”1 

X 

A 

= -0.4686  ft/sec^/rad 

zft 

= -7.269  ft/sec^/rad 

Mft 

= -0.901 6 sec  ^ 

Table  25.  State  Variable  Assignments  for  C-135B  Simulation. 


State 

Variable 

Model  Variable 

Units 

u \ 

9 1 

longitudinal  gust  disturbance 

ft/ sec 

d ) 

w 

g 

vertical  gust  disturbance 

ft/sec 

qg 

pitch  gust  disturbance 

rad/sec 

U 

forward  speed  (body  axis) 

ft/sec 

z 

vertical  displacement 

ft 

w 

vertical  speed  (body  axis) 

ft/sec 

w 

s 

simulator  vertical  speed  (body  axis) 

ft/sec 

6 

pitch  angle 

rad 

q 

pitch  rate 

rad/sec 

9s 

simulator  pitch  angle 

rad 

9S 

simulator  pitch  rate 

rad/sec 

sccl  1 
( 

semi-circular  canal  dynamics 

scc2  ) 

OTO, 

forward  otolith  dynamics 

oto2 

normal  otolith  dynamics 

| 

flight  director  dynamics 

e2  ) 
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Table  26.  Equations  of  Motion  for  C-135B  Model. 

DISTURBANCE  INPUTS 

u = -0.219u  +1.66w, 

9 9 1 

d = -0.592  w + 1.27w0 
9 2 

w = 0.592  d- 1.1 8 w +2.19wn 
9 9 2 

q = 0.00355  d - 0.00708  w - 1 .40  q + 0.0131  w0 
9 9 ^9  2 

AIRCRAFT  DYNAMICS 

0 = 0.03632  u - 0.09922  w - 0.03632  u + 0.09922  w - 32.16489  - 6. 358q 

9 9 

- 0.4686  6 

e 

1 = 0.025  u + 0.9997  w 

w = 0.2553  u + 0.6274  w - 0.2553  u - 0.6274  w + 0.8068  9 + 233  q - 7.269  6 
9 9 e 

9 = q 

q = -0.0006349  u + 0.0007198  d + 0.00246  w + 0.7167  q + 0.0006349  u 
9 9 9 

- 0.003899  w - 0.00098  9 - q + 0.00266  w2 

SIMULATOR  DYNAMICS 

w = 3 . 33  w - 3 . 33  w 
s s 

Ss  = % 
qs  = 2q  - 2qs 


Table  26.  Equations  of  Motion  for  C-135B  Model  (Continued). 


VESTIBULAR  DYNAMICS 
SCC]  = scc2 

SCC2  = -0.001852  SCC,  - 0.08888  SCC2  + 2q  - 2q$ 

OTO,  = -0.2  OTO,  + 0.1990  0 - 0.0932  q + 0.0932  q 

I I S n ^5 

OT02  = -0.2  OTOj  + 0.0207  w - 0.0207  w$  - 0.00502  ^ - 0.7708  q + 0.7708  q$ 

FLIGHT  DIRECTOR  DYNAMICS 

9^  = -0.005  u - 0.1999  w - 2 eh 

0O  = 0.0667  e,  - 0.0667  0_  + 3.82  0 
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The  primary  displays  used  by  the  pilot  for  longitudinal  control  are  assumed 
to  be  the  ADI  (8,9),  ILS  glideslope  (r,  r),  altimeter  (h,  h),  VSI  (Vq,  V^).  In  addi- 
tion, the  Collins  FD-109  Pitch  Axis  Flight  Director  was  available  (Figure  29).  The 
vestibular  inputs  consisted  of  semicircular  sensations  of  pitch  acceleration  and  otolith 
detection  of  specific  force  along  the  forward  and  normal  axes.  The  pilot's  head  was 
assumed  to  be  7.5  ft  above  and  62  ft  ahead  of  the  aircraft  center  of  gravity.  The  dis- . 
played  variable  equations  are  summarized  in  Table  27;  their  corresponding  cost  weight- 
ings and  thresholds  are  presented  in  Table  28.  Only  one  control  input  was  assumed  for 
the  longitudinal  control  task,  since  it  is  evident  from  the  experimental  data  shown  in 
References  58  and  59  that,  as  would  be  expected,  the  throttles  are  used  only  in  the 
very  last  stages  of  the  approach.  Because  the  longitudinal  control  is  through  the  rather 
large  moment  of  inertia  and  stiff  elevator  column,  we  set  the  weight  on  control  rate 
to  correspond  to  a maximum  value  of  approximately  three  inches  per  second. 

5.2.3  MODEL  RESULTS 

Table  29  compares  the  predicted  RMS  longitudinal  errors  with  the  average 

(10,000'  - 5,000')  stationary  statistics  for  selected  variables  reported  in  Reference  59. 

The  model  results  are  presented  for  three  situations,  each  using  a total  fraction  of 

longitudinal  control  attention  of  0.5.  In  the  first  case,  the  attention  allocation  among 

the  five  instruments  was  optimized,  and  the  pilot's  strategy  was  to  spend  as  much  time 

as  possible  on  the  flight  director.  In  the  second  case,  the  control  attention  was 

equally  distributed  among  the  instruments^  = 0.1).  Finally  the  maximum  pitch  rate 

i 

used  in  the  cost  function  was  increased  by  an  order  of  magnitude  to  represent  the  VMC 
situation.  However,  the  results  are  nearly  identical  for  all  three  cases.  Note  that 
the  RMS  pitch  errors  and  RMS  altitude  errors  agree  very  closely  with  the  experimental 
data.  However,  the  experimental  results  exhibit  considerably  more  control  activity, 
and  consequently  larger  pitch  rate  errors  than  predicted  by  the  model. 
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(a)  The  Logic  Symbology  shows  the  "on"  condition. 

(b)  LGS  indicates  glideslope  capture. 

(c)  LGS  45  indicates  45  seconds  after  LGS. 

2.  The  gain  shcedule  is  a function  of  time  from  middle  marker  capture. 


3.  The  gains  and  time  constants  are: 

K.  = Altitude  Rate  Synchronizer  Gain  = 2.0  sec  ^ 

= Altitude  Rate  Gain  --  0.1  deg/fps 
T^  = A I titude  Rate  Lag  T.C.  = 0.5  sec 
T ~ Pitch  Washout  (Lag)  T.C.  = 15  sec 
Kf  = Glideslope  Deviation  Gain  = 16 

Figure  29.  Collins  FD-109  Pitch  Axis  Flight  Director. 


m 
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Table 28.  Display  Variables  - Weights  and  Thresholds. 


Display 

Element 

Maximum  Value 

Cost  Weight 

Threshold 

ft 

5° 

0.04 

2° 

0 

2.5°/sec 

0.16 

1 °/sec 

V 

1/2  scale 

4 

0.1 

Y 

1/2  scale/7  sec 

196 

0.05 

h 

50  ft 

0.0004 

20  ft 

h 

7 ft/sec 

0.0204 

3 ft/sec 

i VD 

200  ft/m  in 

0.000025 

50  ft/min 

^D 

30  ft/min/sec 

0.0011 

7 ft/min/sec 

9C 

4° 

0.0625 

1° 

*c 

2°/  sec 

0.25 

0.2°/sec 

yscc 

- 

0 

0.57  rad/sec ^ 

yoto^ 

- 

0 

0.015  g's 

yoto2 

- 

0 

0.01  5 g's 

Table  29.  RMS  Errors  for  Low-Visibility  Approach. 


Display 

Model 

Experiment  (Ref.  59) 

Variable 

(Units) 

IMC 

Opt.  Attention 

IMC 

= Attention 

VMC 

= Attention 

CAVU 

400/1 600 

100/1200 

© (deg) 

0.73 

0.79 

0.80 

0.77 

0.94 

0.83 

q (deg/sec) 

0.28 

0.29 

0.29 

0.84 

0.70 

0.66 

Y (%) 

10.65 

10.86 

10.81 

h (ft) 

9.77 

9.96 

9.92 

9.58 

8.94 

9.36 

h (fpm) 

90.29 

93.12 

93.05 

6c  (deg) 

0.79 

0.85 

i 

0.85 

6e(deg) 

0.49 

0.47 

0.47 

1.23 

1.26 

1.26 

Independent  research  performed  in-house  at  AFFDL  by  Pollard  (Reference  61) 
and  by  Harrington  (yet  to  be  published)  yielded  similar  results  to  those  shown  in  Table 
29.  These  two  models  plus  the  model  of  this  report,  in  three  independent  studies, 
showed  the  same  trends  for  pitch  rate,  q,  and  elevator  control,  6^.  if  is  theorized 
that  an  incorrect  gain  on  the  order  of  2.25  was  inadvertently  recorded  when  the  experi- 
mental simulation  was  documented. 
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SECTION  VI 


SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

6.1  SUMMARY  AND  CONCLUSIONS 

This  section  presents  a summary  of  the  principle  accomplishments  achieved 
during  this  effort  and  a review  of  the  most  significant  results,  followed  by  some 
suggestions  for  additional  research  and  experimental  investigations  based  on  the  study 
results. 

A comprehensive  literature  survey  was  conducted  to  determine  the  current 
state-of-the-art  of  pilot  modeling  and  to  identify  data  sources  for  developing  and 
validating  the  pilot  model  with  motion  cues.  Section  2 presents  an  overview  of  the 
literature  survey  and  Appendix  A contains  a complete  compilation  of  all  pertinent 
citations  obtained  during  the  study. 

An  extended  pilot  model  was  formulated  to  predict  closed  loop  pilot  behavior 
in  a multi-variable,  multi-loop  aircraft  control  task  with  multiple  inputs  and  outputs. 
The  optimal  control  model  of  the  human  operator  was  selected  as  the  foundation  for  the 
extended  model.  The  conventional  optimal  control  pilot  model  was  extended  to 
incorporate  effects  of  scanning  and  attention  allocation,  motion  cues  sensed  by  the 
vestibular  system,  and  VMC/lMC  (visual  meteorological  conditions/instrument 
meteorological  conditions)  cues. 

There  are  two  dimensions  to  the  approach  of  developing  a pilot  model 
incorporating  visual  and  motion  cues  (indeed,  for  most  model  developments):  one 
dimension  is  the  number  of  data  sources,  ranging  from  raw  data  in  one  well-controlled 
set  of  experiments  to  several  sources  published  in  the  open  literature;  the  other 
dimension  is  that  of  finely  tuning  the  model  parameters  to  provide  excellent  fits  to  the 
data  versus  minimal,  heuristical ly-pleasing  changes  to  the  existing  IMC  models.  The 
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method  we  selected  was  based  on  pragmatic  considerations  of  available  resources  and 
philosophical  grounds.  We  developed  the  components  of  the  model  on  data  published 
by  one  investigator  and  validated  it  on  data  published  by  another  in  a purely  predictive 
fashion.  There  are  many  difficulties  inherent  in  the  use  of  data  published  in  the 
literature,  but  we  feel  that  the  robustness  of  a model  developed  in  this  fashion  will 
compensate  for  any  resulting  lack  of  perfect  fit.  By  the  same  token  we  elected  to 
incorporate  the  components  of  motion  and  visual  cues  by  modeling  their  physical 
properties  and  by  minimizing  the  changes  and  additions  to  the  rules  and  guidelines 
for  specification  of  the  model  parameters.  This  was  done  because  the  limited  number 
of  published  data  sources  would  not  allow  both  a complete  revision  of  parameter 
selection  rules  and  a validation  of  these  numerous  revisions.  The  important  parameters 
of  the  optimal  control  model  of  the  pilot  are  described,  and  guidelines  are  presented 
for  specifying  these  parameters  to  represent  a particular  simulation  situation.  The 
guidelines  we  have  proposed  for  choosing  the  weighting  coefficients  and  thresholds 
removed  much  of  the  unnecessary  freedom  in  the  model,  and  appear  to  give  excellent 
results.  The  guidelines  are  to  choose  the  weights  on  the  display  variables  inversely 
proportional  to  the  square  of  the  maximum  allowable  deviation,  and  to  choose  the 
maximum  display  rate  by  assuming  the  ratio  between  maximum  rate  and  maximum 
position  is  the  reciprocal  of  an  equivalent  time  constant  of  the  desired  response.  The 
thresholds  on  display  position  are  those  derived  from  psychophysical  or  indifference 
considerations,  and  the  threshold  on  the  display  rate  is  that  for  position  divided  by  the 
equivalent  time  constant  of  the  response. 

Motion  cues  may  be  effectively  incorporated  as  additional  measurement  sources 


from  the  vestibular  system  in  the  optimal  control  model  of  the  pilot.  We  have 
accomplished  this  by  modeling  motion  cues  as  being  sensed  by  an  equivalent  set  of 
vestibular  sensors  fixed  in  body  axes  whose  output  is  proportional  to  the  afferent 


firing  rate,  with  the  dynamics  of  the  vestii  jlar  system  incorporated  within  the  controlled 
element  dynamics.  These  vestibular  signals  are  treated  as  additional  measurements  and 
are  not  included  in  the  cost  functional.  Based  on  the  data  used  in  this  study,  we 
tentatively  assign  a noise/signal  ratio  of  -18  dB  to  the  vestibular  measurements.  The 
sensitivity  of  the  model  predictions  to  these  results  should  be  verified  by  a more 
exhaustive  investigation  before  this  value  becomes  an  accepted  constant  in  the  use  of 
the  model . 

VMC  cues  seem  to  be  most  appropriately  modeled  by  different  behavioral 
objectives  based  on  the  data  used  in  this  investigation.  In  the  attitude  control  task, 
we  ascertained  that  the  maximum  roll  rate  used  in  the  weighting  coefficients  for  the 
cost  functional  can  be  increased  by  a factor  of  10  to  describe  the  behavior  of  a 
pilot  using  VMC  cues.  No  data  were  available  for  investigating  the  effects  of  VMC 
cues  in  a position  control  task.  We  suspect  that  these  cues  may  be  accounted  for  by 
increasing  the  maximum  attitude  angle  (threfore  decreasing  its  weight  in  quadratic 
cost  functional).  An  alternative  is  to  consider  an  increase  in  the  maximum  linear 
velocity,  or  a decrease  in  penalty  on  the  velocity  terminal  functional . 

During  the  model  development  phase  three  sets  of  experimental  data  were 
utilized  to  determine  appropriate  values  for  the  model  parameters  to  incorporate  motion 
cues  and  VMC  cues.  The  noise/signal  ratio  for  the  vestibular  measurements  was 
selected  from  using  the  data  base  compiled  by  Shirley  (Reference  17).  Five  controlled 
plants  of  Shirley's  were  analyzed  for  motion  only  (no  visual)  inputs.  It  was  found  that 
three  of  the  five  controlled  elements  gave  a consistent  noise/signal  ratio  of  approximately 
-18  dB  for  agreement  between  the  model  and  the  experimental  data.  The  single 

integrator  (1/s)  gave  higher  RMS  predictions  using  this  value  while  the  neutrally  stable 

2 

plant  10/(s  * 10)  gave  lower  RMS  predictions.  These  results  might  be  explained  by 
the  effect  of  different  levels  of  motor  noise  induced  by  motion  effects. 


The  optimal  control  model  of  the  pilot  was  applied  to  the  VTOL  hovering 
task  described  by  Ringland,  Stapleford  and  Magdaleno  (Reference  23).  These  data 
permitted  the  investigation  of  fixed-base  experiments  in  comparison  with  moving-base 
with  angular  motion  only,  and  moving-base  experiments  containing  both  linear  and 
angular  cues.  Data  were  used  for  simultaneous  longitudinal,  lateral  and  vertical 
tracking  tasks,  with  experimental  results  for  both  "good"  and  "bad1  dynamics 
corresponding  to  the  degree  of  aerodynamic  damping.  The  extended  pilot  model  was 
used  to  predict  RMS  tracking  performance  for  each  task  for  different  levels  of  control 
workload.  It  was  found  that  motion  cues  reduced  the  RMS  error  as  expected,  and 
that  bad  dynamics  generally  exhibit  more  sensitivity  to  changes  in  workload  than 
good  dynamics.  The  optimum  allocation  of  control  attention  between  position  display 
and  attitude  display  was  determined.  Bad  dynamics  were  found  to  require  more  attention 
to  the  attitude  display  to  derive  the  necessary  information  for  stabilization,  whereas 
with  good  dynamics  the  attention  can  be  devoted  to  the  ultimate  position  control 
task.  The  results  of  the  total  task  performance  predictions  showed  excellent  agreement 
between  the  model  predictions  and  the  observed  range  of  subject  means,  both  on  a 
relative  basis  and  on  an  absolute  basis.  The  optimal  attention  allocations  were 
compared  with  the  measured  dwell  fraction  averaged  over  all  motion  conditions.  The 
fraction  of  attention  derived  from  the  model  was  found  to  vary  in  the  same  manner 
as  the  observed  dwell  fraction,  although  there  is  a difference  in  the  attitude  and 
position  displays.  A closer  correspondence  between  the  fraction  of  attention  and 
dwell  fraction  might  be  obtained  by  artificially  imposing  a minimum  constraint 
on  the  fraction  of  attention  for  the  attitude  instrument  to  account  for  pilot  scanning 
habits. 

The  experimental  data  compiled  by  Newell  and  Smith  (Reference  25twere 
used  to  examine  the  effects  of  VMC/IMC  cues  in  the  model  predictions.  For  this 
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roll-control  tracking  task,  the  effect  of  a wide  view  contact  analog  display  (an 
approximation  to  VMC  conditions)  seems  to  be  well  represented  by  an  increase  in  the 
maximum  roll  rate,  rather  than  a decrease  in  roll  rate  thresholds.  Comparison  of  the 
IMC  and  VMC  predictions  with  the  experimental  data  obtained  in  flight  shows  good 
agreement  in  the  mid  to  high  frequency  range.  In  the  low  frequency  range  the  forward 
loop  transfer  function  behaves  as  a gain  instead  of  an  integrator  as  implied  by  the 
experimental  data.  This  low  frequency  behavior  of  the  describing  function  apparently 
cannot  be  resolved  by  simply  adjusting  visual/vestibular  weightings,  but  appears  to  be 
a characteristic  of  the  model  as  it  is  formulated.  (An  alternative  explanation  Is  the 
rational  approximation  to  the  shelf  spectrum  input.  The  optimal  control  model  of  the 
pilot  appears  to  be  quite  sensitive  to  the  input  spectrum  shape,  even  for  the  same  RMS 
level.)  In  spite  of  this  discrepancy  the  agreement  between  the  model  and  the  analytical 
predictions  is  quite  good  ever  the  frequency  range  of  most  interest. 

With  the  parameters  determined  from  the  model  development  phase,  the 
model  was  applied  to  two  problems  for  which  experimental  data  exist  with  motion  cues. 
The  first  experiment  was  conducted  by  Stapleford,  Peters  and  Alex  (Reference  15)  and 
consisted  of  a two  degree-of-freedom  roll-control  task  in  the  NASA  six  degree-of- 
freedom  VTOL  simulator.  The  second  set  of  experiments  consisted  of  data  recorded 
during  a simulated  C-135B  approach  to  landing,  in  which  both  VMC  and  IMC  cues 
were  available,  but  in  which  the  predominant  cues  were  obtained  from  the  instruments. 
In  both  cases  the  model  provided  a very  reasonable  prediction  of  the  experimentally 
observed  results.  The  closed-loop  describing  function  both  fixed-base  and  moving- 
base  agreed  very  well  with  the  experimental  results,  with  the  exception  of  the  above- 
mentioned  fall  off  in  gain  at  low  frequencies.  The  RMS  error  predictions  for  the 
C-135B  approach  also  matched  the  experimental  results  for  pitch  and  altitude  very 


A digital  computer  program  (Progiom  PIREP)  has  been  developed  to  implement 
the  extended  optimal  control  model  of  the  pilot.  The  user's  guide  for  program  PIREP  has 


been  written  os  a separate  volume  for  the  convenience  of  the  user.  This  program  has 
been  installed  and  is  operating  on  the  Wright-Patterson  digital  computer  facility.  This 
is  a flexible  and  useful  program  that  can  predict  closed  loop  pilot  behavior  in  a multi- 
variable,  multi-loop  aircraft  task  with  multiple  inputs  and  outputs.  Program  PIREP  can 
be  a useful  tool  to  improve  simulation  analysis  and  planning,  increase  simulation 
efficiency,  and  lower  simulator  operating  costs.  Furthermore,  aircraft  display  and 
control  system  design,  which  depends  heavily  on  pilot  models,  can  benefit  from  this 
program . 

6.2  RECOMMENDATIONS  FOR  ADDITIONAL  RESEARCH 

The  main  thrust  of  the  recommendations  may  be  summarized  by  the  necessity 
for  validating  and  improving  *he  optimal  control  model  of  the  pilot  under  a wider 
variety  of  motion  and  visual  cue  conditions.  The  specific  needs  are  as  follows: 

1 . More  experimental  data  are  required  to  confirm  the  noise/signal 
latio  of  -18  dB  which  we  have  tentatively  established  using  the 
data  in  this  study. 

2.  Experimental  data  from  VTOL  position  control  under  VMC  conditions 
are  required  to  determine  the  effects  of  VMC  cues  in  this  task,  and 
whether  the  change  in  pilot  performance  can  be  represented  by  a 
decrease  in  the  quadratic  functional  rates. 

3.  A careful  study  on  the  behavior  of  the  optimal  control  model  to  more 
fully  understand  its  low  frequency  predictions  and  the  sensitivity  to 
the  input  spectrum  seems  necessary.  It  appears  as  though  this 
behavior  is  a characteristic  of  the  model  as  it  is  currently  used, 

and  is  not  a result  of  the  incorporation  of  the  motion  cues  used  here. 

4.  The  investigation  of  simpler  representations  of  the  vestibular 
sensor  signals  should  be  undertaken  to  reduce  the  complexity  and 
computation  requirements  of  the  model. 

I I 
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5. 


Perhaps  the  most  extensiw  effort  should  be  expended  in  the 
development  of  quantitative  models  of  subjective  simulator  fidelity. 
It  is  apparent  from  the  many  data  bases  we  examined  in  the  course 
of  the  study  that  the  control  performance  and  describing  function 
are  not  the  most  sensitive  indicators  of  pilot  interaction  with  the 
simulation.  Although  the  with/without  motion  cue  effects  are 
reflected  in  the  describing  function,  the  motion  system  fidelity  as 
represented  by  the  various  washout  circuits  may  lead  to  very  strong 
pilot  opinion  regarding  the  simulation  without  changing  his  control 
performance  or  describing  functions.  Ideally,  such  a simulation 
evaluation  scheme  should  be  able  to  predict  when  the  pilots  feel 
that  the  simulation  is  "good  ',  when  it  feels  "different",  and  when 
the  simulation  is  unnatural,  disorienting,  or  leads  to  debilatory 
response . 
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APPENDIX  A 


LITERATURE  SURVEY  ON  PILOT  MODELING  FOR  MANNED  SIMULATION 

In  accordance  with  the  contract  requirements  a literature  search  was  con- 
ducted. This  included  four  computer-aided  literature  retrieval  searches  in  the  area  of 
pilot  modeling  for  manned  simulation.  A Defense  Documentation  Center  (DDC)  bibliog- 
raphy of  reports  was  searched  for  a ten-year  period  (1965-75)  and  in  the  NASA  Scientific 
and  Technical  Information  Facility  retrieval  of  information  listed  in  STAR  was  similarly 
searched  for  both  limited  distribution  and  open  reports.  The  descriptive  terms  used  for 
both  of  these  searches  are  listed  below: 

Aircraft  piloting  tasks 

Control  theory  pilot  models 

Human  pilot  dynamics 

Manned  aircraft  simulations 

Manned  simulation  of  aerospace  vehicles 

Manual  vehicular  control 

Mathematical  models  of  human  pilot  behavior 
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Pilot  models 
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Although  the  STAR  list  of  open  literature  revealed  little  that  was  not  already  known  to  us, 
a number  of  government  and  industrial  reports  and  theses  not  appearing  in  open  literature 
journals  were  uncovered  by  the  DDC  search. 

An  interactive  computer  retrieval  search  of  the  Compendex  listing  of  Engineer- 
ing Index  citations  was  performed  on  the  MIT  NASIC  system  using  as  key  words:  human 
engineering,  pilots,  simulators,  and  man-machine  systems.  Only  ten  citations  of 
interest  were  brought  forth  from  the  open  literature.  However,  using  the  MIT  NASIC 
system  on  the  NTIS  government  report  data  base  a large  number  of  new  reports  which 
seem  appropriate  to  us  were  revealed.  For  the  NTIS  data  base,  the  following  descriptors 
were  used:  pilot  performance,  simulators,  motion  stability,  vestibular  tests,  and  display 


devices.  We  printed  all  information  on  vestibular  tests  as  well  as  items  having  pilot 
performance  and  any  of  the  other  descriptors  listed.  The  four  literature  searches  were 
then  culled  for  the  titles  and  abstracts  deemed  most  relevant  to  this  contract,  and  were 
divided  into  the  following  main  subject  headings: 
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Pilot  models 
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APPENDIX  B 

OPTIMAL  CONTROL  MODEL  OF  THE  PILOT 


The  overall  closed-loop  pi  lot/vehicle/display  system  is  shown  in  Figure  B-l . 


The  equalization  network  represents  the  means  by  which  the  pilot  attempts  to  optimize 
his  control  strategy  for  a given  situation. 
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Figure  B-l  . Optimal  Control  Model  of  Human  Operator. 


B.l  BASIC  FORMULATION 

To  simplify  the  discussion  the  important  elements  of  the  optimal-control 
model  will  be  reviewed  for  a single  system  control  input  (model  output).  Note,  how- 
ever, that  this  distinction  between  single-loop  and  multi-loop  tasks  is  not  restrictive, 
since  both  are  subsumed  in  this  formulation.  The  multi-control  input  (multi-axis)  case 
is  a straightforward  extension  of  these  results  and  will  be  described  later. 

It  is  assumed  that  the  vehicle  dynamics,  which  may  include  actuator,  sensor, 
and  display  dynamics,  are  represented  by  the  linear  time-invariant  equations  of  motion: 
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x (t)  = Ax(t)  + bu(t)  Ew(t) 


(B-l) 


where  the  n-vector  x(t)  represents  the  vehicle  state,  u(t)  is  the  pilot's  control  input  to 
the  system,  and  w(t)  represents  random  external  disturbances  (e.g.,  turbulence,  wind 
gusts).  Shaping  filters  states  are  incorporated  in  x(t)  to  provide  correlation.  Without 
loss  of  generality,  w(t)  is  a vector  of  independent  zero-mean  normalized  Gaussian 
white  noises. 

Several  system  outputs 

y(t)  - Cx(t)  + Du(t)  (B-2) 

may  be  of  concern,  and  it  is  assumed  that  they  are  presented  continuously  to  the  pilot 
via  some  visual  display  or  instrument  panel.  The  observations  of  y(t)  represent  the 
information  set  upon  which  the  pilot  bases  his  control  actions.  One  may  assume  that 
if  a quantity  y.  is  displayed  explicitly  to  the  pilot,  he  can  extract  the  rate  of  change 
of  that  quantity,  y. . Although  this  has  been  found  to  be  a good  approximation  in 
large,  high-resolution  laboratory  displays,  it  may  be  more  realistic  in  complex  control 
and  monitoring  tasks  to  have  y consist  of  explicitly  displayed  variables  only. 

In  the  model  of  Figure  B-l  an  "equivalent"  observation  noise  vector  is 
added  to  y(t).  A single  noise  v (t)  is  associated  with  each  displayed  output  y.(t). 

Studies  of  controller  remnant  (Reference  B-l)  have  suggested  that  over  a wide  range  of 
tracking  tasks,  the  in'|ected  noises  vy  (t)  are  sufficiently  wide-band  so  as  to  be  approxi- 
mated by  independent  Gaussian  white-noise  processes  with  auto-covariances 

E CV(f>vy/a) } = Vy  • 6 (t-a),  i = 1,  2,  ...,  r 
or 

E tVy(0  Vy’(CT)  } = dia9  (Vy)  * 6 (t_cr)  (B_3) 


1.50 


This  has  been  found  to  be  an  excellent  assi  nption. 


In  general,  a numerical  determination  of  the  noise  covariance  matrix  will 
depend,  among  other  things,  on  the  relevant  features  (quality,  type,  and  form)  of  the 
display  and  on  where  the  human  is  fixating,  i.e.,  whether  he  is  viewing  a display 
directly  (foveally)  or  peripherally. 

The  pilot  perceives 

yp(0  = y(t  - t)  + vy(t  - t ) 

= Cx(t  - t)  + Du(t  - t)  + v (t  - t)  (B— 4) 

a delayed  noisy  replica  of  the  system  outputs.  y^(t)  is  "processed"  by  the  pilot  (through 
some  equalization  network  and  motor  dynamics)  to  yield  a "commanded”  control  input 


ucw- 


A motor  noise  u (t),  which  could  represent  random  errors  in  executing  the 


intended  control  movements,  or  the  fact  that  the  human  does  not  have  perfect  knowledge 
of  the  system  input  u(t),  is  added  to  u (t).  Thus, 


u(f)  = uc(t)  + um(t) 


um(t)  is  assumed  to  be  a (wide-band)  first-order  noise  process  generated  by 


u (t)  + yu  (f)  - u (t) 
mw  ' mw  m'  ' 


where 

u (t)  is  a Gaussian  white  noise  level  with  autocovariance  V . 
m'  ' m 

It  is  assumed  that  the  control  task  is  adequately  reflected  in  the  pilot's 


choice  of  a feedback  control  uc(  •)  which,  in  the  steady  state,  minimizes  the  cost 
functional 


(B— 7) 


J(u)  = E jZjVi  +9Uc 

conditioned  on  the  perceived  information  y^(  • ). 

The  cost  functional  weightings  q.  > 0,  and  g > 0 in  Equation  B-7  may  be 
either  objective  (specified  by  the  experimenter  or  designer)  or  subjective  (adopted  by 
the  human  in  performing  and  relating  to  the  task).  Clearly,  the  selection  of  any 
subjective  cost  weightings  is  a nontrivial  matter  and  is  tantamount  to  mathematically 
quantifying  the  pilot's  control  objectives.  In  some  cases  this  can  be  accomplished  a 
priori,  since  the  actual  weightings  may  often  correspond  to  the  objective  weightings. 
Furthermore,  since  the  pilot  will  adapt  a strategy  to  obtain  "good"  performance,  and 
since  a control  system  designer  would  be  able  to  choose  weighting  to  give  "good" 
performance  in  a linear-system  quadratic-cost  design  procedure,  the  choise  of  weights 
is  relatively  straightforward.  Thus,  the  pilot's  control  characteristics  are  determined 
by  the  solution  of  a well-defined  optimal  linear  regulator  problem  with  time  delay  and 
observation  noise.  It  can  be  shown  (Reference  B-2)  that  the  control  uc(t)  that  minimizes 
J(u),  conditioned  on  the  observations  y^(  • ),  is  generated  by  the  linear  (separable) 
feedback  law 

+ (B-8) 


where  x (t)  is  the  best  estimate  of  the  system  state  x(t)  based  on  the  observed  data 

y (a)  a < t,  and  u (t)  is  the  best  estimate  of  u (t).  Note  that  u (t)  = u (t). 

'p'  " ' m m c c 

The  time  constant  and  the  optimal  gains  L*  are  given  by 


TN  ~ Xn+1 


-1 


Li*  TN  Xi  ' 1 n 


(B-9) 


where 


X = [\j,  Xp/  . ..,  Xn+j  ] is  obtained  from 


l = b„'  K0/g 


(B-10) 


Kg  is  the  unique  positive  definite  solution  of  the  n I-  1 dimensional  Riccati  equation 
Ag’  Kg  + Kg/Xg  + Cg'QCg  - Kgbgbg'Kg/g  = 0 


(B— 1 1 ) 


with  Qg  = diag  [q^ , ^ • • • / qfl;  bQ  = col  [0,  0,  . . . , 0,  1 ] and 


A0 


A [ b 


0 j 0 


Cg  = [C  ! D] 


(B— 1 2) 


The  gain  L * on  u may  be  determined  from 
° m m ' 


Lm*  = TN  I (^1  " A0  + bg^  )”  bg  1 ' (Cg1  + Kgb  )/g  (B-13) 


where 


b ~ col  [b,  0] . 


Two  physically  reasonable  assumptions  are 


Y « T 


-1 


N 


(B-14) 


which  means  that  the  bandwidths  of  uc(t)  and  um(t)  are  approximately  equal,  and 


l * u (t)<<L*x(t)  or  u (t)  — 0 
mm''  ' ' m 


With  these  assumptions.  Equations  (B-6)  and  (B-8)  may  be  added  to  give 


1 Vt 


(B-15) 


TNu(t)  + u(t)  = m(t)  + vm(t) 
m(t)  = -L*x(t) 

The  state  estimate  x(t)  is  generated  from  y ( • ) by  the  cascade  combination 


of  a Kalman  filter  and  a least-mean-squared  predictor,  both  of  which  are  linear  dynamic 
elements.  The  Kalman  filter  generates  a least-mean-squared  estimate  x(t  - t)  of  the 
delayed  state  x(t  - t)  by 

x(t--r)  - A1x(t-T)  + zc0'vy“1  lyp(0-co-(t-T)] 

+ bQTN'1  m(t-  t)  (B-l 6) 

where  x(t)  = col  (x(t),  u(t)  ] . 2 is  the  error  covariance  matrix  and  is  the  unique 
positive  definite  solution  of  the  variance  equation 

0 = AjE  + IA,'  + W - £C0'  Vy-1  CqE  (B— 1 7) 
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The  predictor  generates  the  best  estimate  x(t)  = col  jx(t),  u(t)]  from  the 
Kalman  filter  output  p(t)  = x (t  - t)  according  to 

x(t)  = ?(t)  + eA  1T  |p(t)  - § (t  - t) 

. 5(0  = Aj  5(0  f boTN  1 (B-18) 


1 


Finally,  using  the  techniques  of  Reference  B-3,  it  is  possible  to  obtain  a 
closed-form  expression  for  the  covariance  of  x(t). 


E { x (t)x'(t)  ] = eAlTr  eAlT  + j eAlaWeAl°da 

0 


oo 


+■  f eA  eAl  £ C'V  _1CnV  Al'TeA'ada 
J 0 v U e 

6 


(B— 1 9) 


where  A = A q - b^T  . Thus, 


e r*2  (o } 


Xjj/  i 1/2/  • • • / n 


E [yf  (t) } = (C0XC0')n,  i = 1,  2,  r 


E i>2  (t) } 


= X 


n+I , n+l 


(B-20) 


B.2  MULTI-AXIS  TASKS 


The  description  of  the  optimal  control  model  of  the  pilot  above  is  concerned 
with  single-output  (single  system  input)  models  appropriate  to  multi-group  tracking 
tasks.  The  extension  to  more  than  one  pilot  output,  i.e.,  multi-access  tracking  tasks, 
is  a straightforward  extension  of  the  above  results  (Reference  B-4).  Specifically,  the 
engine  state  Equation  (B-l)  is  modified  to  reflect  multiple  inputs.  The  measurement 
Equation  (B-2)  is  not  modified.  The  cost  function  (Equation  (B-7))  is  augmented  to 
include  control  late  ratings  on  each  of  the  pilot  outputs.  The  control  vector  minimizes 
the  quadratic  functions.  This  is  analogous  to  Equation  (B-8): 

TN  - ^ - = Lx  (t)  (B-21) 
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where  T^  is  a matrix  of  coefficients,  u is  u vector  control  variable,  L is  a matrix  of 
feedback  gains  derived  from  the  deterministic  solution  of  the  matrix  Riccati  equation 
(Reference  B-4)  and  x is  the  output  of  the  Kalman  filter/predictor  in  the  pilot  model. 
When  the  statistics  of  the  noise  are  stationary,  and  when  the  vehicle  dynamics  are 
stationary.  Equation  (B— 2 1 ) reduces  to  the  matrix  of  transfer  functions. 

u(s)  - |Tns  + 1 L x (s)  (B-22) 


B.3  ATTENTION  ALLOCATION  USING  THE  OPTIMAL  CONTROL  MODEL 

The  basic  approach  that  we  follow  in  applying  the  optimal  control  pilot  model 
to  predict  attention  allocation  among  a set  of  display  indicators  is  to  optimize  a quadratic 
cost  functional  with  respect  to  pilot  attentional  constraints.  In  order  for  the  entire  scheme 
to  be  computationally  attractive,  the  following  must  be  accomplished  (Reference  B-5): 

• Relate  attentional  model  parameters  (f.)  to  pilot  model  parameters. 

• Obtain  an  expression  for  J*  = minimum  J that  shows  explicitly  how  the 
f.  affect  the  various  cost  functional  terms. 

dj* 

• Obtain  an  expression  for  the  gradient  terms  that  will  be  needed  in 

subsequent  optimization  algorithms.  ^i 

• Develop  an  algorithm  to  minimize  J*  with  respect  to  the  fj,  subject  to 
total  workload  constraints  on  f.. 

i 

In  this  appendix  the  discussion  is  directed  toward  a pilot  control  task,  although  the 
concepts  can  be  applied  to  a pilot  monitoring  situation  as  well. 

In  the  optimal  control  model,  a fractional  allocation  of  attention  f.  to  the 

informational  variable  y.  modifies  the  "observation"  noise  V (t)  associated  with  that 

^i 

variable  (Reference  B-5).  Thus,  the  noise  covariance  associated  with  y.  is 


lr)6 


N « 2 
v = - or 

i 


V\  f. 


I 


^here 


a.  = o./N(o.) 


(B-23) 

(B-24) 


o.  = RMS  value  of  y. 

i ' l 


N(o.)  = describing  function  gain  of  threshold 


- erfc 


a. 

i 


(B— 25) 


The  noise/signal  ratio  p?  is  the  "full  attention"  noise  ratio,  which  typically  is  O.Olrr  or 


-20dB.  Note  that  Equation  (B-l)  represents  an  implicit  relationship  for  the  actual  noise 


variance  V since  ct?  is  itself  a function  of  V . For  a given  p?/f;,  the  requisite  V is 

y j • y 11  y 


solved  for  via  an  iterative  procedure.  The  quantity  p°/f.  ~ P>  is  called  the  modified 
noise/signal  ratio. 

The  method  that  is  used  to  determine  how  a pilot  allocates  attention  among  the 
various  y-  is  to  minimize  the  optimal  control  cost  with  respect  to  fj  subject  to  constraints. 
This  step  will  re'"  fire  an  iterative  process  to  arrive  at  the  optimal  h,  where  successive 


iterates  f?  result  in  lower  values  for  the  cost  functional.  The  cost  functional  that  is  used 

i 


to  determine  the  fj  is 


J*  - min  J(u) 


(B-26) 


where  J(u)  is  the  basic  cost  functional  in  the  optimal  control  pilot  model. 

In  lieu  of  attempting  to  minimize  fhe  entire  expression  for  J*  numerically  with 
respect  to  f. , it  is  more  efficient  to  isolate  those  terms  in  J*  that  are  affected  by  f. . 


Since  changes  in  f.  are  reflected  as  changes  in  observation  noise  V , we  first  obtain  an 


expression  for  J*  that  shows  the  Vy  dependence.  In  Reference  B-5  it  is  shown  that  the 
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only  term  of  interest  is  called  the  "scanning  cost": 


Jn  - part  of  J * dependent  on  V 


= tr  {L  E L'  } 
1 e e 


(B-27) 


where 


Lg  - "Equivalent"  gains 


Q y2  [T^1  L : 0)  eA 


(B-28) 


and  Z satisfies  the  variance  equation 


0 = AE  + ZA'  +W-EC1  V Cl 

Y 


(B— 29) 


Minimizing  Jq  (and  hence  J*)  with  respect  to  f.  represents  a difficult  nonlinear 
optimization  problem.  The  difficulty  is  two-fold.  First,  f.  affects  in  an  implicit 
manner,  second  affects  Jq  through  the  Riccati  solution  X . As  it  is  unlikely  that  a 
closed-form  solution  for  the  optimum  f.  can  be  found,  the  optimization  process  will  be 
carried  out  numerically  via  some  form  of  gradient  algorithm. 

In  order  that  the  numerical  process  be  reasonably  efficient,  it  is  desired  to 

^ J £ J 

obtain  closed-form  expressions  for  the  gradients  ■ __Q—  or  - - -Q  . Thus,  the  time- 

df;  cWy. 

consuming  process  of  numerically  evaluating  these  derivatives  can  be  avoided. 

We  wish  to  obtain  the  gradient  vector 


(B-30) 


where  is  the  gradient  vector  of  J.  with  respect  to  V _ , i = 1,  2,  ...m,  and  f 


is  a "transformation"  matrix  where 
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(B— 31 ) 


The  gradient  vector  — 9 can  be  evaluated  using  a technique  for  derivatives 
dV 

y 

of  trace  functionals  (Reference  B-5): 


-2.  = diag  (G1  MG) 


(B-32) 


where  the  term  G = FC*  is  the  Kalman  filter  optimal  gain  matrix,  and  M is  easily 
computed  by  solving  the  linear  (Lyapunov)  equation 


A'M  + MA  + L'  L =0 

e e 


(B-33) 


The  elements  of  the  matrix  T must  be  computed.  Since  V = -9_  a . , 

^i  f.  ' 


we  obtain 


dV  . V . i Vy.  V 

r..  = _X  » -Z1  s - 6 ;j  - — sryi 

''  81.  f.  fi  i.  f' 


where  is  a sensitivity  coefficient,  defined  as 
$x  _ dx/x 


If  the  terms  S^'  happen  to  be  « 1,  then  I is  approximately  a diagonal  matrix 


(B-34) 


(B-35) 


r = r..  = 

ii 


(B-36) 


J V) 


More  precise  algorithms  for  detvi mining  the  matrix  r,  including  threshold 
effects  and  additive  as  well  as  multiplicative  observatiort  noise,  are  contained  in 
Reference  B-5. 


SJ 


0 


The  above  algorithm  is  used  to  obtain  g,  = — — 

f df 

gradient  vector.  However,  the  attentional  allocations  f. 


strained  by 


which  is  the  unconstrained 
are  not  free  but  are  con- 


m 

y,  f.  = f*  - total  attention  and  f.  > 0 
M ' 

The  constraint  (Equation  (B-37)  describes  a portion  of  a hyperplane, 
< c , f > = f * 


(B-37) 


(B-38) 


with  c — col  ll,  1,...,!].  Thus,  in  order  to  determine  the  feasible  direction  for 
cost  reduction,  it  is  necessary  to  find  the  projection  of  g^.  on  the  hyperplane  (Equation 
(B-38)).  This  is  given  by 


p __  <gf'  c> 

9f  gf  " c (B-39) 

< c,  c> 


or. 


(B-40) 


of  the  vector.  The  angle  between  g^-  and  g^  is 


Or/  9f  > 
cos  9 = 

Hgf  II  • Hgfll 


(B— 4 1 ) 
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In  order  to  develop  a projected  gradient  optimization  scheme,  assume  we 

are  at  iteration  n,  with  attention  vector  fn . A small  change  Af  such  that  f.  + Af 
still  satisfies  the  constraints  (Equation  (B-3 7))  will  cause  a corresponding  small  change 


in  Jq  . Thus,  at  iteration  n. 


JnH  = I?  + 
J0 


j”  +<gP,  Afn> 


(B-42) 


If  Afn  is  selected  as 


Afn  = “ e JQ  9f/M  9f  II  e<  1 


(B-43) 


then. 


,n+l  \ .n 

J0  -d  -o  J0 


(B-44) 


Thus,  each  successive  iteration  will  result  in  a lower  cost  (to  first  order  terms  only) 
of  100e% . 


Equation  (B_44)  serves  as  the  basis  for  a gradient  optimization  scheme.  We 


set 


e - P cos  9 


(B— 45) 


and  pick  P < 1/2  and  sufficiently  small  such  that  fn  ^ = fn  + Af°  satisfies  f"'  ^ > 0 and 


also  Jq'  ^ < Jq  . Convergence  occurs  when  Jq  ^ is  arbitrarily  close  to  Jq  . If  Jq 
Jq  , a smaller  step  is  taken  by  reducing  P.  Note  that  since  Af  is  in  the  direction  of 
gP,  the  resulting  fn  must  necessarily  continue  to  satisfy  the  constraints  imposed  on  f 
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